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o 

O \ To the memory of Auberon Waugh 

^ ■ (1939 2001), a light in the darkness. 

Abstract 



Observing the Universe, astronomers have concluded that the motion of stars can not 
be accounted for unless one assumes that most of the mass in the Universe is carried 
on by a "dark matter", so far impervious to all attempts at being detected. There 
is now a similar concept of "dark energy". I shall discuss a different subject, "dark 
equations" . These have never indicated that they influence anything or even exist, but 
if one supposes that they do exist, one can systematically discover them and study 
their properties, some of which turn out to be strange and others mysterious. These 
equations are similar in spirit to what one gets when linearizing a given system, or 
studies how an external linear wave interacts with a particular solution of a given 
r^ , system. We define and study linear extensions of dynamical systems in general, and 

^) ' integrable and Hamiltonian systems in particular. Systems discussed include the KdV 

r^ . and mKdV equations and the associated Miura maps, the Burgers hierarchy and the 

associated Hopf-Colc transformations, long wave equations, the Benney hierarchy, and 
the KP hierarchy. 
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1 Introduction 

Politics is for the moment. 

An equation is for eternity. 

Einstein 

Given a system (of equations, say) deserving our attention, we may well inquire whether 
this system is a part of a larger whole still interesting in one way or another. 

In other words, we are talking about extensions, a fruitful mathematical concept. Two 
typical examples are: 1) given a field JF, consider all possible field extensions J^' D J-\ and 
2) given a vector field X on a manifold M, consider (say, vector) bundles ir : E ^ M and 
lifts X ^ X^^* of X from M into E such that 7r*X = X^^V*. If we want to hft up not 
just an individual vector field X but the whole Lie algebra T>{M) of vector fields on M, 
we are looking at connections on bundles with the base M. 

The 2 example is close to what we are after. We start with the following concrete 
problem: 

Problem 1.1. Consider the 2"'^ and 3'''^ KdV equations: 

{X) nt= (3n2 + u^^)^, (1.2) 

(y) Ut = (lOu^ + hul + IQuUxx + Uxxxx)^ ■ (1-3) 

The flows X and Y commute. We seek linear extensions of these flows, X'^^* and Y'^^^, of 
the form 

Ut = (3n^ + Uxx)^, (l-4a) 

ft = Ai{v), (X-*), (1.4b) 

Ut = {Wu^ + 5ul + lOuUxx + Uxxxx)^ , (l-5a) 

vt = A2{f), (y^^*) (1.5b) 

which still commute; here 99 is a vector, and Ai and A2 are matrix linear differential 
operators whose coefficients depend only upon u, not upon ip. 
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We shall discuss the above Problem |1.1| later on, in Section 6. Let us now take a general 
view on this problem. 

It can be looked at from at least four different points of view: 
{A) Given two commuting dynamical systems X and Y: 



(X) ut = X-, 
(Y) ut = Y-, 



find/classify all their linear extensions X°^* and y^^* that continue to commute: 

/yextN 



Ut 
Ut 



X^, ipt = Ai{ip), 



(1.6) 
(1.7) 



(1.8a) 
(1.8b) 



{B) Given an infinite commuting hierarchy of dynamical systems {X^}, find its linear 
extensions {X^^*} that are still commuting. 



(xr) ut = X-, ipt = An{^); 



(1.9) 



{(T) Given the whole Lie algebra of dynamical systems (evolution fields) D^^ = {X}, 
find its linear extensions each of which is isomorphic to -D*^^ as a Lie algebra. In other 
words, if Z = [X, y], for any X,Y,Z € D^^ , then the linear extensions X 
satisfy 



cxt T^cxt '7ext 



yext r x^ext T^extl 

This equality can be further rewritten as 

^[^•^] = X {A^) - Y (^^) - [A^'^A^] , 
where 

(X^^*) ut = X-, ^t = A''iip), 

is the general form of a particular linear extension determined by a map X 
When X and Y commute, the characteristic equation (1.11) becomes 

x(^^)-y(>^^) = [^^,^^], 

which can be thought of as a generalized zero-curvature equation; 
(D) Given a Lie algebra of Ham,iltonian vector fields 



A 



X 



(1.10) 
(1.11) 
(1.12) 

(1.13) 



(1.14) 



find its linear Ham,iltonian extensions. This amounts to finding Hamiltonian matrices 
j^ext^ of the form 

U if 



B 



ext 



B 



6t 



(*) 



(1.15) 



where h is linear in (p. 
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For example, our original Problem LI, with the two KdV flows ( |l.2| ) and (|l.3| ), can be 
looked at as the particular case 

B = d (1.16) 

or 

B = ud + du + ^d^ (1.17) 

of the alternative (D). 

We shall meet every one of the above mentioned 4 cases {A)-(D) later on. At the 
moment let us notice that all these 4 cases are related but neither one covers apriori any 
other. Also, in cases {A)~{C), the solutions {A = A^} can be direct- summed, like (vector) 
bundles over a fixed base M; it's not clear if this observation applies to the Hamiltonian 
case (D), and it's very likely that in general it doesn't. 

The plan of the paper is as follows. In the next three Sections we discuss the general 
case ((7). Section 5 addresses the case (B) for Lax equations, both continuous and discrete. 
In Section 6 we investigate scalar extensions of the first two nontrivial KdV flows ( |1.2D 
and (|l.3|) and find all such extensions, 8 in total. In Section 7 we discuss solutions of the 
extended KdV equations, and in Section 8 we look at Miura maps for these equations. In 
Section 9 we examine extensions given by differential operators of order < 1, for the KdV 
and mKdV hierarchies, and find the associated extensions of the classical Miura map. In 
Section 10 we analyze the Burgers equation and the Hopf-Cole transformation, and then 
generalize all the extensions obtained for the whole Burgers hierarchy. In Section 11 we 
consider the zero-dispersion/zero-dissipation limit of the Burgers hierarchy, and find some 
of its extensions which do not have a dispersive/dissipative origin. Section 12 is devoted 
to the Dispersive Water Waves hierarchy and its quasiclassical limit, while in Section 13 
we study the Benney hierarchy. The Benney hierarchy is the quasiclassical limit of the KP 
hierarchy, and we examine the latter's extensions, two of them, in Section 14. The second 
of these extensions is of a Lie-algebraic nature, and we develop the relevant mathematical 
notion — generalized derivations of Lie algebras — in Section 15. In Section 16 we analyze 
Lie-algebraic meaning of scalar extensions provided by operators of order zero, and arrive 
at nonlinear/relative generalizations of the notion of a derivation of a Lie algebra. In 
Section 17 we arrive at the conclusion that the linear extensions of integrable systems 
are l^*-order approximations to nonlinear extensoins. In Section 18 we find all scalar 
extensions associated with Lax equations with scalar Lax operators of order > 4: it turns 
out that there exist no such extensions in adition to the Lax- type ones and their duals. 
In Section 19 we prove that (nonlinear) Hamiltonian extensions of some general type can 
be composed. 

2 Tangent and cotangent bundles, the functional case 

Let -E be a manifold and q = {qi) be a local coordinate system on E. Let TE be the 
tangent bundle of E, with the induced local coordinates {q,v) = {qi,Vi). Every vector 
field X on E generates a (local) 1-parameter group of diffeomorphisms {At} of E; this 
group lifts uniquely into a group {At} of diffeomorphisms of TE and this latter group's 
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generating vector field X on TE covers the vector field X on E. In local coordinates q, v, 
the motion equations for X and X are, respecively: 



(X) qi = X,{q), (2.1) 

dX, 
dqj 



{X) q. = X,{q), v, = J2^v,. (2.2) 



By construction, the map X i-^ X is a homomorphism of the Lie algebra T){E) into the 
Lie algebra V{TE). This is our first example of the case {C), (iriC|)-(|rT^). 

Turning from Mechanics to Field Theory, we should consider instead of a manifold E — 
a bundle it : E ^ M, instead of vector fields on E — evolution vector fields on the infinite 
jet bundle tTqc : J°°tt -^ M, etc. Rather than rush into the geometric wilderness, let us 
restrict ourselves to the purely algebraic setting; this way we lose topological aspects but 
gain generality by allowing in discrete degrees of freedom and even noncommutativity. 
(Basic facts about differential-difference algebras and related notions can be found in [11] 
and [14].) 

Thus, now the components Xj = Xi{q} = X{qi) of an evolution vector field X are 



functions of {q, x, various x-derivatives and shifts of q}. The lifted vector field X (2.2) 
becomes 

J~) V 

(X) q.=X,{q}, v, = Y^—L(v^), (2.3) 

j -^ 

D _ ( D 



where D = jj- = \ jj- j is the Frechet derivative operator: 



^Df d 

3 ■' 



{f{q + ev)), V/. (2.4) 

e=0 



The map X t-^ X ( p. 3D is an example of (generalized) linearizations [13]. That this 
map is a Lie algebra homomorphism is proven in [13] and [14, Ch. 12] for commutative 
systems and in [16, App. 2] for noncommutative ones. (These references also cover the 
Hamiltonian case {D), but not in the form ( 1.15| ): a different form appears there.) This is 



our 2"^^ example of the general case ((7), and it obviously subsumes the form (2.2). 

To get more examples, we of course could consider any natural bundle over E, not just 
the tangent bundle TE; we won't learn much that way, except proving that the case {(T) 
has an infinite number of solutions. But the case of cotangent bundle T*E will prove to 
be instructive in the next Section, so let's look at this bundle now. 

Let p,q = {p^,qi) be local coordinates on T*E, and let p = pdq = YlP^^qi be the 

i 

universal 1-form on T*E. For any vector field X = '^Xi-^ on E, the expression X\p = 
^p^Xi is well-defined, and the Hamiltonian vector field Xh on T*E with the Hamiltonian 

i 

H = X\p has the following motion equations: 

^^ ^ ■ ^^ ^ -i ^H sr^dXj , , , 
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The map X i-^ ^xjp is well-known to be a Lie algebra honiomorphism; this is our 3'' 
example of the case ( C) . 

Let's now turn to the PDE version of the above cotangent picture. Again, we shall 
stick here — as elsewhere in this paper — to algebraic language. (Geometric version of 
what follows in this Section can be found in [7, Ch. 3].) So, if X is an evolution vector 
field with the motion equations 



(X) qi = Xi{q}, 



(2.6) 



then its cotangent lift is a Hamiltonian vector field Xh with the Hamiltonian H = X\ p 






(2.7) 



The map X i-^ ^x\p is known to be a Lie algebra homomorphism. We thus get a 4* 
example of the case ( C) . 



3 Duality involution 



Suppose we have an extension X i-^ A satisfying the equation (1.11) — for all X in the 
case (C), or for some X in the cases {A) or {B). 
Set 



where f stands for "adjoint". 



(3.1) 



Proposition 3.2. If the family {-4 } satisfies the equation (1.11) then so does the family 

in 



Proof. Applying the operator — | to the equality ( 1.11 ) we get 



A'',A^ 



Thus, extensions come in pairs — if we agree to double-count self-dual extensions. An 
example of the latter is the trivial extension: 



iX 



A^ =0, V X. 



(3.3) 



We have seen in Section 2 two classes of general extensions: the tangent ( |2.3D and the 
cotangent one (|2.7|) . 



Proposition 3.4. The tangent extension ^2.9^ and the cotangent extension (2.1) are dual 
to each other (in the sense of Definition ^.A ))- 
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Proof. For the tangent map X \-^ A (^.31), the matrix A = [A ).. has the matrix 
elements 

U^).. = :^. (3.5) 

^ ^*J Dqj ^ ' 
Hence, by formula ( p.lD , 



Therefore, the corresponding extension has the form 



Comparing this formula and formula (|2.7| ), we see that we need to check the identity 
which follows directly from the definition of the variational derivatives. I 

4 Morphisms of extensions 

Suppose we are given two different evolution vector fields X and Y over two different rings 
Cu and Cy. For example, the field X defines the KdV equation (|1.2|) and the field Y 
defines the modified mKdV equation 

(y) t;i = (-2t;3+T;,,)^. (4.1) 

Let <1> : Cm ^ Cj; be a homomorphism compatible with the vector field X and Y: 

^X = y$. (4.2) 

For example, 

$(ii) = -t;2 + eua;, e = ±1, (4.3) 

in one such (the original Miura map) for the KdV ( |1.2D and mKdV ( [4.1] ) fields. 

The question is, how do such homomorphisms interact with extensions? Clearly, unless 
the map $ is invertible — and thus an isomorphism — only some vector fields on Cu can 
be pushed forward by $; hence, we are in the realm of alternatives (A) and {B) — and 
possibly (D) — but not (C). 

Now, if A provides an extension of X, $ {-^ ) provides an extension of Y: 

(X^^*) ut = X-, ipt = A''i^), (4.4) 

(y-*) vt = Y-, ^t = $ (^^) (^), (4.5) 
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with the extended homomorphism ^^^^ : Cu,ip — > C^^^ actmg as 

^^^\u) = $(u), $^^*(c^) = V. (4.6) 

This mechanism, while general, is however not the only one possible. An interesting 
situation appears when <l>-related fields X and Y are extended tangently, by linearization 
(see Section 2): 

(X-*) u, = X-, ^,= {^), (4.7) 

^ Du 

(y-t) y^ = Y-, ^,= (V;). (4.8) 

Dv 

The general theory [13, 14] then guarantees that the extended fields X^^^ ( ^.7D and y^^* 
(|4.8| ) are related by the extended homomorhpisms <1>'^^* : Cu,ip —^ Cy^^, acting as 

<&-*(«) = cl>(^), q>'^-\^) = £^(^^). (4.9) 

Vv 

(The general theory also guarantees that the map ^^^^ (^]^) is Hamiltonian provided the 
map $ : Cu — > Ct, is. However, this is not applicable to the alternative (D) because our 
form ( |1.15| ) is not of the type furnished by the general linearization theory.) 
For example, the KdV-mKdV case (O), (O), (O) yields: 



(X<=^*) ut = (3n2 + u^^)^ , ift = {Guif + (f^^)^, (4.10a) 

{Y^^') vt = {-2v^ + v^J)^ , Vt = (-6t;V + V'xx),. , (4.10b) 

$<=^t(u) = -t;2 + ev^, $"^*((^) = -2vip + eV'x, e = ±1. (4.10c) 



Of course, ( [4. 5] ) and ( |4.9| ) are not the only kinds of extended homomorphisms possible. 

It is interesting to observe that the duality involution in general destroys extended 
homomorphisms — unless they are of timid type ( [4.5[ ). The tangent point of view, in 
this sense, is vastly richer than the cotangent one — an inversion of the natural order of 
Physics. 

The importance of extended homomorphisms can not be overestimated, for they are 
the most effective organizing principle in studying and classifying extensions. (See an 
example of how this works in Section 8.) 

5 Extended Lax equations 

Lax equations are equations of the form 

Lt = [A,L], (5.1) 

where L and A are operators — in one sense or another. In the differential case [26], L is 
a matrix ordinary differential (or pseudodifferential) operator of the form 

TV 

L = Y^Uid\ d = d/dx, (5.2) 
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satisfying two conditions: 

(i) UN is a constant invertible diagonalizable matrix; 

(ii) un-1 G Ini ad{u]\f). 

(For example, when L is scalar, un = 1 and njv_i = 0.) The general theory [26] then 
furnishes all possible A^s entering formula ( ^.1| ) as 

A = P+, P e Z^{L), (5.3) 

when Z'^(L) is the (abelian center of the) centralizer of L, 

k 

P = Y^p^d^, (5.4) 

with pfc being a constant matrix belonging to the (center of the) centralizer Z'^(ui\f) of mat; 
the + and — notation employed in formula (|5.3[) means: 



[Y^asdA =J]a,9^; (j^as^M = J] a.S^ (5.5) 

The main result from the general theory [26] we shall make use of below is: 



Denote by Xp the evolution derivation (5.1) with A = P+: 



XpiL) = [P+,L], PeZ%L). (5.6) 

If Q is another element of Z^{L), 

Xq{L) = [Q+,L], (5.7) 

then 

Xp{Q) = [P+,Q] (5.8) 

and 

[Xp,Xq]=0. (5.9) 

With minor changes, similar results apply to discrete Lax equations [8, 11]. Here L is 
a scalar operator 

AT 

L = ^UiA\ UN = l; (5.10) 

— oo 
A is the (dual) shift acting on functions of n G Z: 

{A'{f)){n) = f{n + s), n,seZ; (5.11) 

P = L"^, m G N; 

(^a.AM =5]a,A^ (j^a.AM =^asA'^. (5.12) 



All other formulae (^.1|)-(|0|) remain unchanged. 
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Now, let P run over the generators of Z'^{L), and let Xp be, for each P, the corre- 
sponding evolution derivation (|5.6| ). We are, therefore, in the alternative {B) of Section 1. 
We define the following linear extension of the derivations Xp: 

Lt = [P+,L], (5.13a) 

^t = P+iv). (5.13b) 



Proposition 5.14. All the extended flows X'p^ (5.1c.) commute between them,selves. 



Proof. Let P,Q a Z^{L). By formula ( 1.13 ), we need to check that 

Xp(Q+) - Xq(P+) - [P+, Q+] = 0. (5.15) 

This is one of the Wilson's formulae [26], and it follows from formula (|5.8| ) and the following 
calculation: 

Xp(Q+) - Xq(P+) - [P+, Q+] = [P+, Q]+ - [Q+, P]+ - [P+, Q+] 
= [P+,Q+ + Q_]+ + [P+ + P_,Q+]+-[P+,Q+] 
= [P+,Q_]+ + [P+,Q+]+ + [P_,Q+]+ 
= [P+ + P_, Q+ + Q-]+ = [P, QU = 0, 

since P and Q commute, being elements of an abelian center of the centralizer Z'^{L). ■ 

Remark 5.16. Let K = 1 + ■ ■ ■ be the dressing operator for L, conjugating L into its 
highest term: 

L = KuNd^K~^, (5.17a) 

L = KA^K''^, (5.17b) 



for the differential and difference cases respectively. The Lax motion equations ( |5.6| ) result 
from the Wilson motion equations [27] in the dressing language: 

dp{K) = -P_K. (5.18) 

Then the extended system ( |5.13| ) results from the extended system 

Kt = -P-K, (5.19a) 

^t = P+iv). (5.19b) 



Obviously, all the extended systems ( 5.19 ) still commute in the larger dressing space. 

Remark 5.20. Set 

L = -L\ P = -P\ (5.21) 

and apply the duality involution to the system ( ^.131) . We get: 

Lt = [P+,L], (5.22a) 

ipt = P+{ip), (5.22b) 

which is a system of the same type as ( [5.13 ), though not necessarily the same system: 
( 5.13| ) and (^.22|) are dual to each other. 
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Example 5.23. Let 

L = d'^ + u (5.24) 

be the KdV Lax operator. Then P = (-^^") , n G Z+, and therefore 

L'^ = L, P'f = -P ^ L = -L, P = P, (5.25) 

and all the extended KdV flows ( 5.13| ) are self-dual. 



Remark 5.26. What does the duality involution do to the extended dressing system 
( 5.19| )? Denote by Lq the highest term of L entering the equations ( 5.17| ), so that 

L = KLqK-^. (5.27) 

Taking adjoint, we get 

L = kLoK-^, (5.28) 

where 

K = K''^^. (5.29) 

Therefore, by formula ( |5.19a ). 



Kt = -K^-^ (k^\ i^t-i = _i^t-i (-K^p}\ ii't-i = p^K = -P_K. (5.30) 

Thus, the dual form of ( 5.1S| ) is: 



Kt = -P-K, (5.31a) 

ipt = P+i^). (5.31b) 

Remark 5.32. The Lax equation ( ^.13a| ), 

Lt = [P+,L], (5.33) 

results from the compatibility analysis of the Lax eigenvalue problem 

L{ip) = Xip, (5.34a) 

supplemented by the time evolution of the eigenfunction ip: 

n = P+{ip). (5.34b) 



The Lax system (|5.34 ): 



L{ip) = Xif, ift = P+{p) (5.35) 

is not the same as our "extended Lax system" ( 5.13| ): 



Lt = [P+,L], ^t = P+{^). (5.36) 



Our system (5.36) is more general, for it does not assume that (p is an eigenfunction of 
anything; in other words, the constraint 

{^'^^L{ip) = const} (5.37) 

is compatible with our extended system ( ^.36|) but it is not implied by that system. 
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6 The first two KdV flows 



Let us now tackle the Problem 1.1. We need to make some assumptions to make this 
problem manageable. We make three such assumptions: 



(a) Lp is scalar, 

(b) The extension operators Ai ( |1.4b ) and A2 ( |1.5b| ) are polynomial in 
and x-derivatives of u; 

(c) The extended systems preserve homogeneuity of the original ones. 



(6.1) 



The homogeneuity count is based on the ranks 

rkfu'^'A =i + 2, rk{d) = l, (6.2) 

rk{X) = 3, rk{Y) = 5, (6.3) 

where u^^' = ^-~. Thus, we look for the KdV extensions of the form: 

ut = d (3n^ + Uxx) (6.4a) 

ift = {aua, + (3ud + jd^) (c^), (6.4b) 

ut = d (Wu^ + 5ul + Wuu^x + u^^A , (6.5a) 

ift = {{auxxx + buux) + [cuxx + dv?) d + euxd"^ + fud^ + gd^){<p), (6.5b) 

where a, (5, 7, o, 6, c, d, e, /, g are unknown constants. 

We assume, in this Section, that Ai and A2 are of maximal order so that 7 7^ and 
hence 5 / 0. The case {7 = 0=^g = / = e = 0}is treated in Section 9. 

Equating the like-terms in the equation 

X{A2)-Y{Ai) = [Ai,A2l (6.6) 

we arrive at the following relations: 

(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 
(6.7e) 
(6.7f) 



/ = 


37 


e = 


i^+«- 


d = 


3 7^ V 37 


c = 


37 V 3 


b = 


3 72 V 3; 


a = 


5<7 2 + 7-1 
37 3 
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d 
b-- 



30/? 
12-/?' 
60a 

12-/3' 



(6.8a) 
(6.8b) 



/? = 



9^ 



7 



1) (7-1 + 2) + (5 - 1) 



0, 



(5- l)a = (7-1)0, 

3{fa-al3 + b) = a(c+ 10), 

a{f - 10) + a(6 - /?) + 6(1 - 7) = 0, 

/?(/ - 10) + c(6 - /3) + 2d(l - 7) - 376 + 3a f = 0, 

18c - 20^ - 6^7 - 97/3 + (e + c)p + 2ea = 0, 

(e — c — /)a + a/3 = 0. 



(6.9) 

(6.10a) 
(6.10b) 
(6.10c) 

(e.iod) 

(6.10e) 

(e.iof) 



We break the analysis of the system (|6.7|) - (|6.10| ) into four steps, as follows. 
Step 1: a = /? = 0. Then a = b = c = d = e = f = 0, and we get the decomposed 
systems 



ut = d (3n^ + u^x) , ift = ^d^{ip), 



(6.11) 
(6.12) 



Obviously, such decomposed extensions are available for any system of vector fields any- 
where. For the KdV hierarchy, we then have 



{X"^^) ut = X^, (/?j=7„a2" ^((^), nGN, 7„ = const, 

and all these are self-dual extensions. 
Step 2: P = 0, a ^ 0: no solutions. 
Step 3: a = 0, P ^ 0: there are three possible solutions: 

Q = 0, /3 = 6, 7 = 1: ^t = {Qud + d^){ip), 

and this is dual to the linearized case ( |4.10a| ); 

a = 0, /? = 3, 7 = 1; / = e = c = 5, d = 10, 5 = 1, a 
ut = d {3u^ + Uxx) , (pt = {3ud + d^) (if), 

ut = d (lOu^ + 5ul + lOnn^.^. + n^''^' 

^t = {{5uxx + lOn^) d + 5n^a2 + 5ud^ + d^){ip); 

a = 0, /3 = -6, 7 = -2; / = e = -20, d = c = -10, 

g = -4, a = 6 = : 

ut = d (3n^ + Uj^j.) , (pt = - {Qud + 25^) (99), 



(6.13) 



(6.14) 

(6.15) 
(6.16) 

(6.17a) 
(6.17b) 

(6.18) 
(6.19) 
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ut = d ( lOn + 5u^ + lOuUxx + u 



(4) 

ift = -(10 {uxx + u'^)d + 20uxd'^ + 20ud^ + 4d^){ip). 

Step 4-' oifi 7^ 0. There are four possible solutions: 

a = 3, /3 = 6, 7 = 4; a = 15, 6 = d = 30, c = 50, 
e = 60, / = 40, 5 = 16 : 



(6.20a) 
(6.20b) 



(6.21) 
(6.22) 



Ut = d (3m^ + Ui:^^) , (ft = (3m(? + 3du + 45^) (c^o), 

nt = 9 (Wu^ + 5ul + Wuuxx + u^'^A , (6.23a) 

ipt = ((15u^.^.^. + 30uu^) + {bOuxx + 30^2) d + 60u^.5^ + 40m9^ + I6d^){ip), (6.23b) 
and this is the self-dual Lax form ( |5.24| ), ( ^.25| ); 

a = /3 = 6, 7 = 1; a = / = 10, 6 = 60, 
c = e = 20, d = 30, g = 1 : 

this is the linearized case (|4.10a ), dual to ( |6.14| )~( |6.17| ): 

Ut = d {3u^ + Uxx) , (pt = {d{6u + d'^)){ip), 
ut = d (Wu^ + 5ul + Wuuxx + u^'^A , 

(^4 = {{lOuxxx + QOuux) + {20ux + SOu^) d + 20uxd^ + Wud^ + d^){ip) 
= {d {{SOu"^ + lOuxx) + Wuxd + lOud"^ + d'^}){ip); 



a = (3 



7 = 1; a 



20, 



d 



10, / = 5, 



9 



1 



Ut = d {3u'^ + Uxx) , ipt = {d{3u + d'^)){ip), 

ut = d (Wu^ + 5ul + WuUxx + u^'^^' 
ft = {{5uxxx + 20uux) + {Wuxx + lOu^) d + lOu^O^ + 5ud^ + d^){f) 
= {d { {5uxx + Wu^) + 5uxd + bud"^ + d'^}){ip), 

and this is dual to the case (|6.15| )-( |6.171 ); 

a = 13 = -6, 7 = -2; a = d = -10, 6 = / = -20, 

c = -30, e = -40, g = -4: 

ut = d (3m^ + Uxx) , ^t = {-d (6u + 29^))(y?), 

ut = d (Wu^ + hu\ + Wiuuxx + u^'^A , 

ip^ = -{(Vduxxx + 20to^.) + [3^Uxx + lOu^) d + 40m^9^ + 20Ma^ + 49^) (c^) 
= -(9 {10 [uxx + u^) + 20M^a + 20^9^ + 49^})(v3), 



(6.24) 

(6.25) 
(6.26) 



(6.27) 
(6.28) 
(6.29) 

(6.30a) 
(6.30b) 



(6.31) 
(6.32) 

(6.33a) 
(6.33b) 



and this is dual to the case ( 6.1S| )-( |6.2C| ). 

Thus, we have found 8 solutions: the decomposed one ( p.ll| ), ( |6.12 ) and the Lax 
one (|6?2lD -( |633|) , both self-dual; the linearized one (|6^ -( |6^26r and its dual ( |6ll) ; the 
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strange one (|6.28D ~ (|6.30D and its dual ( 6.15| )~ (|6.17| ); and the mysterious one ( 6.3l| )~ (|6.33| ) 
and its dual (|6l8|)-(|63o|). 

Which ones among these 8 extensions are applicable to the whole KdV hierarchy and 
not just to the two nontrivial KdV flows? 

The decomposed extension is certainly applicable, and is given by formula ( |6.13| ). The 
Lax extension likewise applies to the whole hierarchy, by formulae ( 5.13| ). Ditto the lin- 



earized extension, by the results of Section 2, and its dual, by the results of Section 3. 

This leaves us with two mutually dual pairs. The mysterious extension (|6.31| )~ (|6.33| ) 
appears at the moment just that, mysterious. (It will be explained in Section 17.) Let's 
look closely at the strange extension ( |6.28|) - (|6.3O0 . Let us compare the strange and the 
linearized extensions: 

ut = d (3n2 + u^^) (KdVa), (6.34a) 

ift = (69n + (9^) (if) (linearized), (6.34b) 

ipt = {3du + d^) (if) (strange), (6.34c) 

ut = d (lOu^ + 5ul + Wuu^^ + u'^'^A (KdVs), (6.35a) 

ipt = {d{ (30u^ + lOu^^) + lOn^a + lOud'^ + d"^} (((/?) (linearized), (6.35b) 

^t = {d{{Wu'^ + 5u:^x) +5u^d + 5ud'^ + d'^}){ip) (strange). (6.35c) 

We see that in each of the two cases, the strange operator ylstrango can be gotten by the 
formula 

^strange — p^ (D.oDj 

from the vector field ut = X^, where the integrated operator D = -^ acts on differential 
(-difference) polynomials in u (without constant terms) via the rule 



Du 



-^ = f dt^{tu}, /{An} = AV{^} V A = const, (6.37) 

r Du Jq Du 

where r > is the total u-degree of a homogeneous polynomial /; for / non- homogeneous, 
the definition ( |6.37 ) is extended by additivity. It is very likely that formulae 



D^X'^ 

ut = X-, ^t = —-JL{^), riGN, (6.38) 

Du 

provide strange extensions for the whole KdV hierarchy. 

In the next two Sections, we shall approach the KdV case from different perspectives. 

7 Solutions of the extended KdV equations 

We start with general linearized systems: 

Proposition 7.1. Suppose X is an evolution vector field invariant with respect to trans- 
lations in x".' if 

Ut = X- (7.2) 
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is the corresponding motion equation then its every solution u = u{x, t) generates solutions 
u{x + Te'^^t) V T, where 

e° = (0, ... ,1,... ,0) (1 atthe a*'' place). (7.3) 

Let 

ut = X-, ipi = ——{^p) (7.4) 

Du 



be the linearization of X. If u is a solution of ( 17. ^J then 

is a solution of ( |7.^ j. 

Proof. Translation invariance assures that if ut = X'^{u} then 

^{u{x + Te",t))=X-{u{x + Te'',t)}, V r. (7.6) 

ot 

Differentiating this with respect to r at r = we get: 

d f du\ DX^ f du\ 

(7.7) 



dt ydx'^J Du V5x" 
Thus, (^ solves {^. m 

This is as far as general linearization results go. Let's look now at the particular case 
of the KdV equations. 

We normalize the KdV flows by requiring the highest u-derivative in X„, u^'^"'~^' , enter 
with coefficient 1: 

(X„) Ut = d{hn), hn = n(2«-2) + . . . ^ ^ ^ pj_ (7,8) 

The differential polynomials hn are in fact variational derivatives of some Hamiltonians Hn'- 

K = ^, (7.9) 

du 

and it is known that 

{d'^ + 2ud + 2du){hn)=d{hn+i), n(^'L+, (7.10) 

/lo = 1/2, /ii = n,... (7.11) 

because the KdV equations are bi-Hamiltonian. 
If we write 
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then formula ( 7.10| ) yields 



AuOnnvT' ^Ux + 2ux9nu"- = On+i{n + l)n"na: 
2(2n + l)^ 



7n+l 



n + 1 
2«-i(2n-l)!! 



nl 



formula ( 7.12| ) tells us that (for n > 1) 



Pn 

(2n + l)pn 



™,2n+3 



4{c + n(n + 1)} 



(7.13) 



n G N. (7.14) 



If we look at stationary solutions of the n KdV equation, of the form 

u = ex \ c = const, (7.15) 



Mn = M-{2n-2) (7.16) 

^ M = -2. (7.17) 

Thus, 

u = c/x^. (7.18) 

To find all possible values of c, set 

pi = c. (7.20) 

Using formula ([7.10|) , we find 

{-2n{2n + l)(2n + 2) - 4c • 2n - 4c} 



2(2n+l), . ^. . ^ 

^ p„+i = ^-— -^{c + n{n + IK 7.21 

n + 2 

^ Pn(c) = ^—. ^-\{{c + ^{l + l)}. (7.22) 

Thus, all stationary solutions of the n^^ KdV equation of the form u = cjx^ have the 
c-values 

c=-^(^ + l), £ = 0,...,n-l. (7.23) 

By Proposition |7.1| , the corresponding linearized flows are satisfied by 

99 = const/x^, n > 1. (7.24) 



By Corollary 7.35 below, the dual to linearized flows are all satisfied by 



93 = const/a; , n > 1. (7.25) 
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Proposition 7.26. Suppose a Hamiltonian system 

ut = d{h), /i = -— , 
du 

has a stationary solution u = f{x) such that 
= 0, 



6H_ 

5u 



u=f{x) 

and the corresponding linearized system 

ut = d{h), (ft = dj^{(p), 
has a stationary solution {u = f{x), (p = u;(x)} such that 

{u^{x)) = 0. 



Dh 



u=f(x) 



Then the dual to linearized system {7.2i): 



Dh 



DuJ 



Ut = d{h), (ft -- 
has a stationary solution {u = f{x), ip = 0,{x)}, where 
n{x) = [ uj{C)dC. 



Proof. By the basic property of variational derivatives [20], 



Su J \ 5u 



so that 



Dh\^ _ Dh 
15^) ~ IJu' 



Hence, for the RHS of the (/9-equation ( [7.31 ) we get 



DhV 
IDu) 



u=f(x) 



d{Q{x)) = ^ 
Du 



u=f{x) 



[by (^] 

(w(x)) ^^-^ 0. 



(7.27) 



(7.28) 



(7.29) 



(7.30) 



(7.31) 



(7.32) 



(7.33) 



(7.34) 



Corollary 7.35. If u = f{x) is a stationary solution of ( 7.27 ) and H is x-independent 



then {u = f{x), if = const f{x)} solves the dual to linearized system ( 7.31 ). 
Proof. By Proposition |7.1|, we can take uj{x) = f'{x). 
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By formula ( 7.23 ) with 
u = -21 x^ 



1, 



(7.36) 



is a stationary solution for each KdV flow j^n with n > 1. Let us see what stationary 
solutions of the form Lp = const x^ look like for the 7 extended KdV2 flows of the preceding 
Section (from the total of 8, with the decomposed extension deleted). 
The linearized extension ( 4.10a ) 

a (6n + 9^) ((/?)= : A^=-3,4; 

The dual to the linearized extension ( |6.14 ) 

(6u + 9^) 5(^) = : iV = -2,5, 

in agreement with Corollary 7.35| ; 
The self-dual Lax extension ( |6.22| ): 

(3iia + 39m + 49^) (99) = : A^ = -l,l,3; 

The strange extension (|6.29D : 



(7.37) 



(7.38) 



(7.39) 



d (3u + 9^) (^) = : N 



-2,3; 



The dual to the strange extension ( |6.16 ): 
(3m + 9^) a((/p) = : iV = -l,4; 

The mysterious extension ( |6.32 ): 

-5 (6n + 25^) (^) = : Af = -2,3; 

The dual to the mysterious extension ( 6.1S| ): 



(e-u + 2^2) 5((^) = 



N 



-1,4. 



(7.40) 



(7.41) 



(7.42) 



(7.43) 



We see that the sets of strange and mysterious exponents are identical, an artifact of 
the 2""^ extended KdV flow where 



x^ext 



■^^ strange 



M- 



(7.44) 



'-myst(¥') 

No similar simple relation exists for the 3*" extended KdV flow; but a more complex one 
may. 



Remark 7.45. Proposition 7.2(: and Corollary 7.35 can be considerably generalized, as 
follows. 
Let 



--'f) 



(7.46) 
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be a Hamiltonian system whose Hamiltonian matrix B is g-independent. Let 

"_D_ /my 



..^B,f), „,.B 



(7.47) 



and 



«.-(f). ^. 



D (5H 



B{p) (7.48) 

be the corresponding (tangent and cotangent) hnearized extension and its dual. 



Proposition 7.49. (i) If 

{q = f, p = g} 



is a solution of the cotangent system l \7.4W then 
{q = f, v = Big)} 



(7.50) 



(7.51) 



is a solution of the tangent system ( 7.4V ; 
(a) The map ^: 

q = q, v = B{p) 

is a Hamiltonian map between the canonical Hamiltonian structure 



(7.52) 



J 



1 
-1 



(7.53) 



)/ the cotangent system ( 7.4i ) and the linearized Hamiltonian structure [14, P- 192] 



B 



-B 
-B 



(7.54) 



of the tangent system (7.4''i)- 



Proof, (i) Apply the operator B to the 2"^^ equation in ( 7.4^ ). This results in the 2 
equation in ( 7.47 ); 

(ii) The Jacobian of the map $ ( [7.52 ) is 



nd 



Therefore, 

Jac ■ J ■ [Jacy 



(7.55) 



1 W 1 
B )\ -1 

-B 
-B 



-B 



B. 



(7.56) 
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On the other hand, the corresponding Haniiltonians are also related by the map ^ ( [7.52 ): 
the linearized Haniiltonian [14, p. 182] is 



(v) ~ —V ■ —— = —B{p) ■ — — ~ —p ■ B' I — — = p ■ X, 



Dq 



6q 



6q 



6qJ 



where 



^^B,f) (.„) 



is the original evolution vector field (7.46|). 



Proposition 7.59. Let 



Ut = d{h), h 



6H_ 

6u ' 



be a Hamiltonian system with the Hamiltonian H that is x -independent. If 
u = f{x,t) 



is a solution of ( \7. 6(\ ) then 
{u = f, p = const /} 
is a solution of the cotangent system 

Ut = d{h), pt = jT-dip). 
Du 

Proof. Since h is j;-independent. 



Ut = 


cJ{h) 


so that 




ft = 


Dh 
Du 



Dh 



Ux) 



Dh 

mi 



diu), 



d{f), 



u=f 



and this is the form of the p-equation in ( [7.63| ) . 

The conclusion ( [7.62 ) applies also to some non-Hamiltonian systems: 
Proposition 7.66. Let 

ut = X^: Uii = X^, i = l,... 



(7.57) 



(7.58) 



(7.60) 



(7.61) 



(7.62) 



(7.63) 



(7.64) 



(7.65) 



(7.67) 



be a dynamical system where each X'^ is a quasipolynomial, i.e., a sum of quasimonomials, 
the latter being products of the terms of the form 



u 



{^) 



±1 



(7.68) 
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in other words, of terms of positive and negative degrees. Assume that no X^ contains 



terms of total degree zero, and consider the following strange extension of ( 7.6V ): 

D^X' 



ut = X^, ipt 



Du 



-i^)-- 






(7.69a) 
(7.69b) 



where 



D'f _ 1 Df 

Duj deg (/) Duj 

for any quasihomogeneous f of total degree deg (/) ^ 0. 
If u = f{x,t) solves ( 7. 6'^) then 



{u = f{x,t), ip = const f{x,t)} 
solves ( 7.6!^ . 
Proof. The statement ( [7.71 ) amounts to the identify 






Dui 



f 



(7.70) 



(7.71) 



(7.72) 



for any quasihomogeneous / of deg (/) 7^ 0. 

Now, 



/=(nM)', O^iGZ 



D^f / - ^x^-i 



DUn 



i', 4"' 



«"-ESm-/- 



(7.73a) 
(7.73b) 



Finally, if / and g are quasihomogeneous, with / deg (/), deg (g), deg (/) + deg (g), and 
formula ( 7.72| ) holds true for both / and g, then for h = fg we find 



D^h 



1 Dh 



1 D{fg) 



Duj deg (h) Duj deg (h) Duj deg (/i) \ D 
deg{g)f—— + deg{f)g 



,Dg Df 



m-i 



deg {h) 



Dui 



Dui 



E ^(^^O = xrrTTTtdeg {g)fg + deg {f)gf] = h. 



Dui 



deg {h) 



Du-, 



Corollary 7.74. // / is quasihomogeneous then 
E;^(%)=d«g(/)/' deg(/)GZ. 



(7.75) 
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Proof. If deg (/) 7^ 0, this is just a reformulation of formula ( 7.72| ). If deg (/) = 0, we 



decompose / as / = gh, with deg (g) deg (h) ^ 0, and then use the derivation property of 
the Prechet derivative operator jj:^- B 

It seems likely that, for every fixed scalar differential Lax operator L, 

n-l 

L = d'' + Y,Ujd''-^-\ (7.76) 

1=1 

all the commuting Lax flows 

Lt = [P+,L], P£Z{L), (7.77) 

still commute after being strangely extended. 

8 Miura maps for extended KdV equations 

The KdV equation 

ut = (3n^ + Uxx)^ (8.1) 

and the mKdV equation 

vt = {-2v^ + vxx)^ (8.2) 

are related by the classical Miura map 

u = -v'^ + evx, e = ±l. (8.3) 

We saw in Section 4 that the linearized equations 

ut = d (3n2 + u^^) , ipt = d{6u + d^) {^) , (8.4) 

vt = d {-2v^ + V,,) , ^Pt = ^ {-Qv^ + 8^) (V) (8.5) 

are related by the linearized Miura map (4.10cD: 



u = —v'^ + evx, if = —2v'iIj + eipx, e = ±1. (8.6) 

Let us see which general extensions of the KdV 

ut = d (3n^ + Uxx) , ^t = {aux + /3ud + 7^^) (c^) (8.7) 

and of the mKdV 

vt = d {-2v^ + Vxx) , (8.8a) 

^t = {{avxx + bvvx + cv^) + [dvx + ev"^) d + fvd'^ + jd^) (ip) (8.8b) 

are related by a generalized Miura map which, from dimensional considerations (assuming 
it's polynomial), must have the form 

u = -v'^ + evx, ip = ei/jx + Ovi; = {ev + ed){ij). (8.9) 
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Here a, (5, 7, a, 6, c, d, e, /, are unknown constants, with 7^0. 
Calculating ipt in two different ways, from (|8.7| ) as 



V^i = {q(-2ot^ + ei;^^) + /? (-t;^ + ef^) 9 + 75=^} {Ov + ea)(V), (8.10a) 

and from (|8.9|) as 



934 = {6* (-277^ + V^^)^ + {Ov + e5)((0f:j;:j; + ftVf^j; + cv^) 

+ (dt;^ + ev"^) d + /vS^ + -fd^)] (ip), (8.10b) 

and equating the expressions ( ^.lOD , we find the following relations: 

/ = 0, (8.11a) 

e = -P, (8.11b) 

d = €{P + 3'y6), (8.11c) 

c = 0, (8.11d) 

h = 2{(3-a)- 876!^ (8.11e) 

a = e(Q-/?); (8.11f) 

0(7^2 _^^ 2) = 0, (8.12a) 

0/3 = 2(/? - a) - 87^^ (8.12b) 

6»(7 - 1) = Q - /?, (8.12c) 



If = 0, the relations ( 8.12| ) collapse into 



a = /? = arbitrary, (8.18a) 

and this case is not interesting, for all it says is that 93 is a conserved density and one can 
introduce the potential ifj such that 

if = ei>^. (8.18b) 

So, suppose 

e^O. (8.14) 

Then (^.12a| ) and ( ^.12c ) can be considered as providing (3 and a in terms of 7, 9, while 



J.12b[ ) is the desired relation on 7 and 9: 



p = 2 + -f9\ = 2 + 7^2^^(7-1), (8.15) 

= -9f3 + 2(/9 - a) - 87^^ = e{-2 - 76*2 + 2(1 - 7) - 876*} 

= -9-f{e^ + 39 + 2) = -9-f{9 + l){9 + 2). (8.16) 

Since ^,7 / 0, we obtain 

9 = -1,-2. (8.17) 
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Thus, 

^ = {ed + ev){i^) = ( ^'j^/^J ) i^P) for ^ = ( 12 ) ' (8-18) 

an interesting and unexpected result. 
Summarizing, we get 



so that formulae ( |8.7[ )-( p^ become 

^t = {3u^ + {2 + 7)ud + jd^){ip), ^ = {€d-v){i^), (8.21a) 

i^t = ((e(l - 7)v.x - (2 + 7)^^«x) + (2e(l - 7)^^^ - (2 + 7)^') d + 79^)(V'), (8.21b) 

<ft = {{^ + 2j)u^ + {2 + Aj)ud + jd^){v), ^ = ie^-2v){^l^), (8.22a) 

^j = ((2e(l - 7)7;^^ - (4 + 87)1;^;^) 

+ (2e(l - ^)v^ - (2 + 47)^2^ d + jd^){^). (8.22b) 

Let us determine which ones of the 7 KdV extensions of Section 6 come out of the 
corresponding mKdV ones. 

From formulae ( |8.19| ) we see that: 

a = <^ 61 = -2, 7 = -2 => /? = -6, (8.23) 

and this is the dual to mysterious extension (|6.18| )- (|6.2[l|) ; 

9 = -l, 7 = 4, a = 3, P = 6, (8.24) 

and this is the self-dual Lax extension ( 3.2l|) -( |6.23|) ; 



e = -2, 7 = 1, a = /3 = 6, (8.25) 

and this is the expected linearized case (|8.4D -( p^ ; 

9 = -l, 7 = 1, a = /3 = 3, (8.26) 

and this is the strange extension (|6.28| )-( |6.3[l| ), with 

n = -rri^)^ V = TtW' ' 8.27a 

Du Uv 

i;t = {-^vv^ -2,v^d + d^)W + ^ -^ ^^ (V^) . (8.27b) 

uv 
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The latter formula shows that the strange extensions of the higher niKdV flows, if they 
exist, are not given in terms of the intergrated operator ■^. 
The Lax case ( ^.24 ) deserves some comment. The formula 

ip = {ed-v)ii;) (8.28) 

is certainly unexpected, for Lax equations and modified Lax equations have essentially the 
same eigenfunctions (see [29]). Now, with the KdV Lax operator 

L = d'^ + u, (8.29) 

ip satisfies the eigenvalue problem 

L{ip) = Xip, (8.30) 



which becomes, in terms of ip ( ^.28 ): 

Lied-v){ip) = X{ed-v){tl;). (8.31) 

Recall [19] that the Miura map (|1|) 

u = -v^ + ev^ (8.32) 

comes out of factorization of the Lax operator d^ + u: 

d'^ + u={d-ev){d + ev). (8.33) 

Hence, the equation (|8.3l| ) becomes 

{d - ev){d + ev)e{d - ev){ip) = Xe{d - et;)(V'), (8.34) 



or 



{d + ev){d-ev){i;) = X'ip, (8.35) 

which can be rewritten as 

L(^) = A^, (8.36) 

where L is the Backlund transform of L: 

L = d^ + u, u = -v^-ev^. (8.37) 

The time evolution of tp ( 8.21] ) is 

^i = (3 {-ev^ - v^)^ + 6 {-ev^ - v^) d + 4d^){^P) 

= {3ud + 3du + Ad^) (ip) = P+ {ijj) , (8.38) 

with 

P = 4L^/2_ (839) 
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Thus, every Lax extension 

Lt = [P+,L], ^t = P+{^), PeZ{L), (8.40) 

comes from the corresponding modified Lax extension 

vt = ..., i^t = P+W (8.41) 

under the extended Miura map 

u = -v'^ + ev, ip= {ed-v){'ip), L = -v'^ - ev^. (8.42) 

Similar conclusion applies to all scalar Lax equations with the Lax operator of order n >2: 

n-2 

L = d'^ + ^Uid\ (8.43) 

with the Miura map being (differential) factorization 

n 

L={d + vi)---{d + Vn), J^^^i = 0, (8.44) 

and with the Backland transformation being effected by a cyclic permutation of the roots 

v/s: 

L = {d + V2)---{d + Vn){d + vi), (8.45) 

^ = {d + vi){ilj), (8.46) 

^t = P+{^), 4^t=P+W. (8.47) 

9 KdV and mKdV extensions of order < 1 

In Section 6, we classified commuting KdV extensions (|6.4| ), (|6.5| ) 

ut = d {3u^ + Uxx) , ipt = {aux + /3ud + ^d^) (if), (9.1) 

ut = d (Wu^ + 5ul + Wuua:x + u'^^A , (9.2a) 

ipt = {{auxxx + huux) + (cuxx + du^) d + euxd"^ + fud^ + gd^){ip), (9.2b) 

The assumption there was that 'jg ^ 0. Let us now consider the case 

7 = 0. (9.3) 

It's immediate that 

g = f = e = 0, (9.4) 

and then 

c = /3, d = /?(/? + 4)/2, (9.5a) 

P{P - 2) = 0. (9.5b) 
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If 

/3 = => a = 0, a = b = c = d = 0, (9.6) 
and we get the trivial extensions. So, let 

p = 2 =^ c = 2, d = 6. (9.7) 
Then 

a = a, b = 6a, a is arbitrary, (9.8) 

and we finally obtain: 

ut = d (3n^ + Uxx) , (9.9a) 

<^t = {au^ + 2ud){ip), (9.9b) 

ut = d (lOu^ + hu\ + lOnn^.^. + n^^^) , (9.10a) 

^t = {a {u^x + 3u^)^ + 2 (n^^ + Sn^) a)((^). (9.10b) 



We are now going to explain formulae (|9.9| ), ( |9.10 ), and then determine similar exten- 
sions for the whole KdV hierarchy. 

Proposition 9.11. Let X and Y be two commuting evolution vector fields, with the mo- 
tion equations 

{X) qt = X-, (9.12a) 

(Y) q, = Y-. (9.12b) 

Consider their scalar extensions, of the form 

(X^^*) qt = X-, ipt={v + V) {^), (9.13a) 

(y^^*) qt = Y-, ^t={w + Wy^), (9.13b) 

where 

a a 

are vector fields, andv, w, V°^ , W" are all functions of{q}. 
Then the extended systems ( 9.13a ) and ( 9.13l\ ) commute iff 



X W] -Y V 



\V,W] : (9.15a) 



X(w) - Y{v) = V{w) - W{v). (9.16) 
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Proof. This is just equation ( |1.13 ): 



X iw + W] -Y iv + V 



v + V,w + W 



written out for the summands of order 1 and 0. 



(9.17) 



Corollary 9.18. (i) When V and W satisfying ( 9.15 ) are fixed, the solution space {v,w) 
of equation ( \9.1(\ ) is a vector space; 

(a) When at least one of V and W is not divergence-free, the dimension of the vector 
space in (i) is >1. 

Proof, (i) is obvious; 

(ii) Set first v = w = Q, a solution of ( p. 16 ). Thus, 



iX 



\Y 






A^ = V, A' = W. 
Applying the duality involution to expressions ( |9.19|) , we find 

A^ = - (^^)t 

Thus, 

V = div ( V 
is a solution of ( |9.16| ). 



dw(V]+V. 



w = div ( 1^ j 



(9.19) 



(9.20) 



(9.21) 



Corollary p.l8| explains the zero-order terms in equations ( 9.9b| ) and ( |9.10b| ) , aux and 
Oi {uxx + 3n^) , in terms of the l**-degree terms 2u and 2 {uxx + 3n^) , respectively. 

Let us now interpret the l^'^-degree terms themselves. Comparing formulae ( |9.9| ) and 
(|9.10| ), we see that the extended KdV fields are of the form: 

[XZ') Ut = d{hn), ipt = iahn-l,x + 2hn^id){ip), (9.22) 

where the differential polynomials hn are those introduced in Section 7. By formulae ( |7.9D , 
( 7.10| ), we can rewrite ( |9.22 ) as 



ud + du + (9^/2 —(fx + -^zd^p 



"Px +2"^^ 







6/6u 
6/6ip 



(H, 



n-l 



(9.23) 



Since the Hamiltonians {Hn} commute in the 2 Hamiltonian structure of the KdV 
hierarchy 



B^^ = 2{ud + du + d^/2) 



(9.24) 



the extended KdV flows ( |9.22D commute as well, provided the matrix entering ( |9.23| ) is 
Hamiltonian. That it is so indeed, can be seen after rewriting one-half of this matrix as 



B 



ud + du (pd — Xdip 
dip — Xipd 



+ 



d^/2 




— 1^ 



(9.25) 
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The linear part of the matrix B corresponds to the semidirect sum Lie algebra Vi k V , 
where Pi is the Lie algebra of vector fields on M^, and V^ is the one-dimensional Pi -module 
of A-forms {function x (dx) }; the constant part of B ( |9.25| ) is a generalized 2-cocycle on 
Pi; thus, S is a Hamiltonian matrix (see [14]). 

It remains to show that the thus constructed KdV extensions ( 9.221) exhaust completely 
the centralizer of Xf^* ( |9.9|) . So, suppose 



(X^^*) ut = dihn), ipt = n{ip) (9.26) 

commutes with Xf^^: 

X2{n) - Xn{2ud) = [2ud,n]. (9.27) 

This is an inhomogeneous equation on TZ, with one solution being given by formula ( p. 22 ). 
The corresponding homogeneous equation is then 

X2{n) = [2ud,1l]- (9.28) 

If 

N 

n = Y.fsd\ fs = fs{u}, (9.29) 

s=0 



then equation ( 9.28 ) becomes 

N N N 

Y,X2{fs)d' = J22ufs,.d' -2j2fs[d'',u]d. (9.30) 

s=0 s=0 s=0 

Picking out the 9^-terms, we get for Jn = f- 

X2U) = 2ud{f) - 2fNu,. (9.31) 

If / is independent upon {u} then / must vanish. If / is not independent upon {u}, let i 
be the highest x-derivative of u entering /. Then 

^2(/) = E ^ (6-^^^ + -''f = ^-t.^-'^'' - ^^-'"'f- (9-32) 

jr=o i=o 



The term 

-u^'+-'^ (9.33) 



^/ -M+^) 



on the LHS of ( |9.32| ) is not matched by anything else. Thus, the only solution of ( p. 31 ) 
is / = 0, and so 7^ = 0. (The same argument works for a linear combination of X!j!^^'s 
(PI).) 
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Remark 9.34. The Hamiltonian form ( |9.23| ) provides an extension of any commuting 
hierarchy with the Hamiltonian form (|9.24| ), not just the KdV hierarchy. We have here, 
thus, an example of the alternative {D) of Section 1. 

Let us now examine the case 7 = (|9.3|) from the point of view of Miura maps studied 



in Section 8. Equation ( ^.16 ) is now automatically satisfied, and we find: 

p = 2, a = 2-e; e = -2, d = 2e, h = 29, a = -€9: (9.35) 

nt = a (3n2 + -u^^) , >^t = {{2-e)u^ + 2ud){^), (9.36) 

vt = d (-2t;3 + v^^) , (9.37a) 

i't = {-e {ev^ - ^') , + 2 {ev^ - v^) 8) (V) (9.37b) 

= {-eu^ + 2ud){'tp), (9.37c) 

u = ev^-v'^, Lp={td + ev){'4j), (9.38) 

is arbitrary. (9.39) 

Thus, the reduced KdV extension ( |9.36| ) comes out of the modified KdV extension ( p. 37]) 
under the Miura map ( tJ.38D . In the old notation (|9.9|), 



a = 2 -6. (9.40) 

We now make a bold leap of faith and declare that the same Miura map ( |9.38| ), when 
applied to the following mKdV„ extension 

vt = d{td + 2v){hn^i), i^t = {-0K-i,x + 2K-id)iij), (9.41) 

produces our KdV„ extension ( |9.22| ): 

ut = {d^ + 2ud + 2du){hn-i), ^t = {{2-e)K-i,^ + 2hn-id){v). (9.42) 

Proof. Call hn~i by h. Since the 2"*^ Hamiltonian structure d^ + 2ud + 2du of the KdV 
equations is related to the Hamiltonian structure —d of the mKdV equations by the Miura 
map u = eVx — v'^, the motion equations 

ut = {d^ + 2ud + 2du){h), h = 6H/du, (9.43) 

and 

vt = -d{6H/6v), (9.44) 

are related by the Miura map 

u = evx-v'^. (9.45) 

But 

6H/6v = Jac\6H/6u) = {ed - 2u)\h) = -{ed + 2v){h) (9.46) 

^ vt = {-d){-€d-2v){h) = d{ed + 2v){h). (9.47) 
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This justifies the Vf-equation ( |9.4l|) . 

Now, 

^t = [{ed + ev){ij)]t = 9vt^ + {ed + 0t;)(Vt) 

= {Ovt + {ed + ev){-eK + 2hd)){ij), (9.48a) 



while ( p. 421) predicts that 

(^j = ((2 - e)K + 2hd){^) = ((2 - 0)/i^ + 2hd){€d + 0v)(V'). (9.48b) 

Thus, we have to check that 

Ovt + {ed + ev){-eK + 2hd) = ((2 - e)h^ + 2hd){ed + Ov), (9.49) 

and this identity is readily verified. I 



Notice that the extended Miura map ( 9.38 ) is no longer Hamiltonian w.r.t. the Hamil- 



tonian structure ( |9.23| ). (This can be most easily seen by letting 9 vanish.) 

Remark 9.50. Our analysis of KdV extensions, at least for the KdV equation itself, is 
now reasonably complete. However, that analysis was made under the assumptions (a)-(c) 



(6.1). These assumptions may be too restrictive. Assumption (a), that (p is scalar, is most 
evidently so. The niKdV hierarchy is associated with a 2 by 2 eigenvalue problem, so a 
pair of eigenfunctions enters there. It seems that in general, for scalar Lax operators of 
order n: 

n-2 

L = d'' + Y,Uid\ (9.51) 

i=0 

the independent building blocks of all possible extensions have the sizes 1,2, .. . ,n. (All 
intermediate sizes are provided by the number of factors of partial factorizations of L.) 
See more on scalar extensions in Section 18. 

Assumption (b), that extension operators A are polynomial in {u} (but not rational, 
say), is also suspect. Recall [29] that at the root of the KdV theory lies the Mobius- 
invariant Ur-KdV equation 

(9.52) 



(9.53) 



*i 


3 

= -^xxx - 2* 


L/*x, 




which is 


connected via 


the potential 


map 


P = 


--^x 






with the equation 






Pt 


= pxxx - - {pi 


/P)x' 




which, in the variable 






w - 


= linp, 







(9.54) 



(9.55) 
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takes the form 

wt = Wxxx - '^wl, (9.56) 

which is the potential form 

V = Wx (9.57) 

of the mKdV equation 

vt = {vxx - 2v^)^ , (9.58) 

which is mapped via the Miura map 

u = evx- v^ (9.59) 

into the KdV equation 

ut = {3u^ + Uxx)^- (9.60) 



The equations ( p. 521 ) and ( |9.54 ) are themselves rational, so the polynomial assumption on 



their extensions would be totally unreasonable; by implication, the polynomial assumption 
(b) for the KdV and mKdV equations may be not as natural as it had appeared. 

Assumption (c), that extensions are homogeneous, seems secure at the moment. We 
shall see later on, in Section 12, that it has to be modified to remain true in spirit, if not 
in form. 

10 The Burgers hierarchy 

The Burgers equation, 

ut + uux = i^Uxx, V = const, (10-1) 

discovered by Bateman [1], can be rescaled into any convenient form desirable. We choose 

ut = d{u^ + Ux). (10.2) 



The famous Hopf-Cole transformation makes this equation (10.2) out of the heat equation 



Vt = Vxx ■■ (10.3) 

u = Vx/v. (10.4) 

We adopt in this Section the following useful language. Let C{f } be the differential 
field generated by an indeterminant v = v^^' and its x-derivatives {v^-^'lj € N}. This 
field has a differential subfield C{u} C C{i;}, with u ( |10.4| ) being v^^' /v. The evolution 
derivation X2 : C{v} -^ C{v} dTo!^ ), 

X2iv) = v^^\ (10.5) 
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restricts onto subfield C{n}: 

X2{u) = 2nn(i) + n^^) = d (u'^ + u'^^A . (10.6) 

One, thus, can introduce the whole Burgers hierarchy [3], by starting with an infinity of 
commuting evolution derivations X„ of Clti}: 

X„(t>) = !>("), neN, (10.7) 

and then restricting these derivations onto the subfield C{n}: 

Xn{u) = d{Qniu)). (10.8) 

Here 

Qn{u) = {d + ur{l) (10.9) 

are differential polynomials in u, satisfying the defining relation 

v^""^ = vQn (v^^^/v) = vQniu) : (10.10) 

If 

V^''^=vQn, Qn = Qn{u), U = V^^^ /v, (10.11) 

then 

Qo = l, Qi = u, Q2 = u^ + u^^\ ... (10.12) 

and 

= v{uQn + d{Qn)) = vQn+1 (10.13) 

^ Qn+1 = (5 + 7X)(Q„), (10.14) 

and formula ( 10.9]) follows. 

We aim to determine all possible extensions of the Burgers hierarchy subject to the 
restrictions ( |6.lD : we are looking at extensions that are scalar, polynomial, and homoge- 
neous. 

It will later on prove very useful to start with the full f -picture first. All extensions of 
( 10.7] ) can of course be immediately written down: 

(X^"*) vt = v^''\ V'i = en^^"\ e„ = const, n € N. (10.15) 

Let us estimate what kind of subfields of C{v,ip} the flows ( [L0.15| ) could be restricted 
onto. If e„ = 1, the flow X^^ can be thought of as a linearization of the flow Xn (10.7). 
Therefore, the restriction C{i;} D C{u}, 

u = d{inv) = v^^^v, (10.16) 
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can be extended by linearization into 

^ = ^^i^l^) = ^{i;/v). (10.17) 

Dv 

To calculate cpt, we introduce temporarily a new variable (^ such that 

ip = d{ip), (10.18) 

V' = vif. (10.19) 

Therefore, 

^(n) =±(1) .("-^)^W = . E (I ) Qn->^^''' (10-20) 

fc=0 ^ ^ fc=0 ^ ^ 

^ ^(")/t; = Qn^ + J2{l) Qn-kV^'-'^- (10.21) 



Hence 



ipt = d{{'4)/v)t) = d 



V V V 



= d I enQn^ + ^n^(l] Qn-W^^'^^ " ^Qn J (10.22) 

^ ^t = d{{en- l)Qn^ + e- Z] ( ^ ) Qn-kV^^~^ J ■ (10.23) 

We see that we must have 

e„ = 1, (10.24) 

and thus, the flows 

(Xr) vt = v^'^\ V'i = V'("), (10.25) 

under the restriction 

u = v'^'^^/v, ip = d{'ip/v), (10.26) 

become 

(X-t) ut = d{Qn), (10.27a) 

^^ = (9T.{l) Qn-kd'-'){v>)- (10.27b) 

Since the extension ( 10.27| ) is in fact linearization, we get: 



^ = E ( : ) Q"-^'-^ (10.28a) 



fc=i 

n 



-C^Qn o v^ / n 



Dit ^-^\ k 

k=Q 



^ = E 1 Qn-fe^'^-Qn. (10.28b) 
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Analysing the above calculation, one can observe that it will still go through if we 

(10.29) 



replace the restriction ( 10.17 ) 

ip = d{i)/v) = V'^^V'^ - ^i^^^V"^^ 

by only "half" of it. Set 

(p = tp^ > /v, 
so that 

ipy ) = vip 

Therefore, 



(10.30) 



^in-k)^ik) ^ ^ 



k=0 



E 

fc=0 



n 



Qn-W 



(k) 



(10.31) 



u—n V / 



t'y^^''^ - QnV ■ 



n 



{-Qn + en ^ ( 1 ) Qn-kd^ 
fc=0 ^ ^ 



fc=0 



>M- 



(10.32) 



This restriction is different from the previous one, ( 10.27b| ). 

If we now go over the derivation of formula (10.32), we may notice that the restriction 
ip = ipy^i Jv ( |10.30| ) can be further generalized, into 



f 



-(A^A 



(10.33) 



where p is arbitrary (formal symbol). We should extend the differential fields C{v} and 
C{v, ■0} by adjoining t;^, but this is a minor matter since 



(7;^)(l) = pvPv^^^/v = yPpU. 

Consequently, 

(Z;P)W = vfQr^ipu). 

Hence, from formula ( |10.33| ) we obtain: 



(10.34) 
(10.35) 



(,P^) W =±[1] {vn^-^^ ^('^ =-'E(1) Qn-kipu)^^'^ (10.36) 

A:=0 ^ ^ fc=0 ^ ^ 



Vt 



fc=0 ^ ^ 
{-pQn{u) + en ^ ( ^ jQn 



-k{pu)d^){^). 



(10.37) 
(10.38) 
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The thus obtained extension ( 10.38| ) yields, for the 2 and 3'' Burgers equations, in 
the notation 



£2 = 7> es = 9, 
respectively: 

ft = { {pilP - l)u'^ + p{j - l)ux) + 2'j pud + 7^^) {(f) , 

ut = d {u^ + 3uux + Uxx) , 

ft = {{p {gp^ - l) n^ + 3p{gp - l)uux + p{g - l)uxx) 
+ 3(7 {p\^ + pux) d + Sgpud^ + gd^){y^). 



(10.39) 



(10.40a) 
(10.40b) 

(10.41a) 
(10.41b) 



Are there other nontrivial subfields of Clv^ip}, besides the two we have used, ( |10.26[) 
and ( |10.33 ), that are invariant under extended Burgers hierarchy? Let us see what can be 
said about the 2"*^ extended Burgers equation 



Vt = Vxx, Ipt = li^xx, 

restricted by the general form 

onto the system 

ut = d {u^ + Ux) , ¥?* = {{aux + uju"^) + Pud + jd'^){ip). 



(10.42) 



(10.43) 



(10.44) 



Here A{v) and B{v) are unspecified functions of v to be determined. A straightforward 
calculation shows that there exist two types of restrictions: ( |10.33| ) for general 7, and, for 

7 = 1, 



Vt = Vx 



-p-1 



Vx + V~ 



'd) (tP) = v-P{d - u){ij) = v^-Pd{i;/v), 



ut = d {u^ + Ux) , 

ipt = {2ux + p{p - iW + 2/Jn(9 + 9^) ((/?). 

The linearization map ( [lO.lT] ) is the case {p = 1} of formulae ( |10.45| )"( |10.47 ). 



(10.45) 
(10.46) 

(10.47a) 
(10.47b) 



To see that the subfield C{u, (f} C C{v, t/j} (|10.45) is invariant under arbitrary flow 



(^r) Vt 



,(«) 



i^t = i^ 



(n) 



let us temporarily introduce a new variable (/? such that 

ip = v(p 



(10.48) 



(10.49a) 
(10.49b) 
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Using formula ( |10.21| ), we obtain 



(1 - p)Qn^ + 



k=l ^ ^ 



Now, 






(10.50a) 
(10.50b) 

(10.51) 



((l-p)Qn + E(^) 



{d+ip- l)u)Qn-kid + {p- l)uf-')i^). (10.52) 



For p = 1, this formula turns into formula ( 10.27b| ), while for n = 2, formula ( |10.47b| ) is 
recovered. 

Our goal now is to determine all commuting extensions of the 2 Burgers equation, 



(I10.44D , and the 3^"^ one: 

Ut = d {u^ + 3uUx + Uxx) , 

(ft = {{auxx + buux + cu '^ + [dux + eu ) d + fud + gd )(</?), 



(10.53a) 
(10.53b) 



where a, uj, /?, 7; a, b, c, d, e, /, g are unknown constants. We shall assume for the time 
being that 7 7^ 0, and defer the nonmaximal case 7 = until later. 

Writing out in long-hand the equality ( 1.13 ) for the flows ( |10.44| ) and ( 10.53 ) and 
equating the like-terms, we find: 



UJ 



PiP - 2) 

47 



27 



472^ ' 



d = f^[2a + /3(l + 7-^)], 
47 



c = g 
b = - 



P 



P 



2-fJ 27' 
3P 3 g 



3 g 



2747^ 27^ 



-1 + 



(/3 + 2)(l + 7-^) 



a = -^ + ^{3a + P)-'-^l'-[2a + P{l+j-% 



(10.54) 
(10.55a) 
(10.55b) 
(10.55c) 

(10.55d) 

(10.55e) 

(10.55f) 



a 



7 



1 



27 



Pa 



,^1-7^3 



(10.56) 
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Gguj — 3jb + da + b + 6a — ga = 0, 

-26 + 6uj + da- 3af + 3/3a = 0, 

2foj - 67c -Pb + 2duj + 46 - Gw - 6a = 0, 

2fuj - (3b + 2ea - 2duj + 6c - 6a; = 0. 

Equation ( 10.56| ) says that either 
7-1 



a 



27 



-/3 



or 



g-l + 



1-723 



72 4 



75 = 



(10.57a) 
(10.57b) 
(10.57c) 
(10.57d) 



(10.58) 



(10.59) 



(or both). Equality ( 10.58 ) is covered by formula ( 10.40b| ). 

So, we are left to tackle the relation ( 10.59| ), and of course the 4 relations ( 10.57 ); 
formulae ( 10.55 ) we treat as defining the constants /, e, d, c, 6, a in terms of a, /3, 7, g. 

Our strategy in handling these opaque equations is to look at them as perturbations 



around the known solution (10.40b). In other words, we set 



9 = 13/2-1 => /3 = 2p7, uj = p{jp-l), 
a = p{-f -I) +7A, 



(10.60a) 
(10.60b) 



so that A is just [a — p{-~f — l)]/7. The values /, e, d, c, 6, a from formula ( 10.41b| ) we 
designate as "old" : 



/ 



old 



^9P, 



e"'' = 3gp\ 
d""" = 3gp, 
c°"^ = P {9P' - 1) 
6°'^ = 3p(5p-l) 
a°id = p(5_l). 



(10.61a) 
(10.61b) 
(10.61c) 
(10.61d) 
(10.61e) 
(10.61f) 



These old variables satisfy the relations ( |10.57] ). The new variables, read off formulae 
(I1055D and (|10.60|) , are: 



a = a^'' + Mi±2) A = Pig 



47 



1) + Mill) A, 

47 



6 = 6°'^ + M2^±l) A = 3p{gp - 1) + M2^±I) A, 



old 



27 
2 



27 



c = c"'" = p {gp^ - I) , 

d = d°'^ + ^gA = 3gp + ^gA, 

e = e"'^ = 3gp\ 
f = r''' = 3gp. 



(10.62a) 

(10.62b) 

(10.62c) 

(10.62d) 

(10.62e) 
(10.62f) 
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Now, equation (|l0.57dD is satisfied identically for any A. Equation ( |l0.57c ) yields: 



= (4 - 27p)M^;^±^A + 2p(7P - l)^gA - 67A 

= 3A I ^ (2 + 7p - 7V') + gp{pi - 1) - 27! = 3A 
= 6A7-1 {g - 7^) . 
Since we are interested in A 7^ 0, 
5 = 7^ 



7 



Substituting this into ( 10.59 ) we find 

7' = 5 = 1- 
With the latter relation granted, formula (|10.57a) yields: 



(10.63) 
(10.64) 
(10.65) 



= (1 - 37)^(/o + 7)A + 3pA7 + ^Ap{j - 1) + ^AA7 + 6 • ^(1 + 7)A - 9A7 

„. .1 — 37 7(1 — 37) 7 — 1 A7 3,^ N „ 

3A <! ^-^p + ^^^-^ + 7P + /5^ + ^ + 2 (1 + 7) - 37 



3A(-7 + ^ 



3A7 ( - + 1 



(10.66) 

^ A = 2. (10.67) 

The difference between equations ( |10.57a| ) and ( 10.57b| ) is linear in A, and is identically 

satisfied for 7^ = 5 = 1 ( pM5| ). 

Thus, the complete list of extensions ( 10.44 ) with nontrivial centralizer ( 10.531) is: 
(|10.40b|); 



7 



a 



P = 2p, uj = p{p-l), 



(10.68) 
(10.69) 



and this is the case ( 10.47 ); and 

7 = _1, a = 2/5 -2, 13 = 2p, u; = -pip-l), 
which is dual to (|10.68D : 

-{2u^ + p{p - 1)^2 + 2pud + dy = (2/3 - 2)u^ - p{p - l)v? + 2pud - d"^ . (10.70) 

The extension ( |10.69| ) is the only one on this list which doesn't come from the restriction 
of an evolution derivation from C{v, if)} onto C{n, (/?}; however, its dual ( 10.68| ) does. This 
suggests that the family ( [10.40b| ) is self-dual; indeed, the whole family ( 10.38| ) is: 

Proposition 10.71. Denote by An = Anip, e^) the operator entering the RHS of equation 



An{p,(-n) = -pQn{u)+en^Yli\ k ) ^"■ 



-k{pu)d\ 



Then 



-An{p,eny = An{-p,{-ir+'en) . 



(10.72) 



(10.73) 
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Proof. We need formula ([12, p. 64]) 

n 

E 



k=0 



Qn-k{pu)d^ = {d + puY 



Then, 

An{p, e„) = -pQn{u) + e„(5 + /ju)", 
and therefore 

-An{p, e„)t = pQ^{u) - en{-d + puT = pQn{u) + {-ir+\n{d - pur. 
Let us now examine the reduced case 7 = 0: 

ut = d [u^ + u^) , ift = {{aux + uju^) + Pud) (if). 
It's immediate that f = g = in (|l0.53b| ), so that 

93j = {(auxx + 6uux + cu ) + (dtii; + eu ) d){ip). 
Equating the 9^-terms in the defining identity (1.13|), we find 



e = d = P, 
/3(/?-l)=0. 

We postpone the case /? = 'till later. So, suppose 

P = l = d = e. 
By Corollary 9.18 and formula ( |9.2l| ), we can take 

a = 0. 
This will force 

a = 0. 
The remaining part of the equation ( 1.13| ) yields: 

b (nX~)^ + 3cu X~ — 2u!uY^ = u [buux + cu ) — [ux + n ) • 2u!uux, 
where 

X^ = d [u^ + Ux) , y = 9 {y? + 3uux + Uxx) ■ 



(10.74) 



(10.75) 



(10.76) 
(10.77) 



(10.78a) 
(10.78b) 



(10.79) 
(10.80) 



(10.81) 



The terms UxUxx in ( |10.81 ) provide 6 = 0; then the nn^-terms force w = 0; and the terms 
u^Uxx then lead to c = 0. Thus, 

a = u; = 0; /3 = 1; a = 6 = c = 0, d = e = l (10.82) 

is the essential solution. It can be generalized for the whole Burgers hierarchy: 
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Proposition 10.83. The extended flows 

{X^^') ut = d{Qn), ^t = {Qn-id){ip), nGN, (10.84) 

all commute. 
Proof. We have to show that 

Xn{Qm-l) — Xm{Qn-l) = Qn-lQm-1 ~ Qm-lQn-V (10.85) 

Now, 

JJu Uu 

[by ( |l0.28b|) , (|ia7^)] . _i r^ \(r^\ nnoR\ 
^ ((a + n)™ ' -Qm-i) (Q„), (10.86a) 

while 

Qn-lQi^Ll = Qn-l{{d + U)- n)(g„_i) = Qn-liQm " uQ^_i). (10.86b) 

Hence, the equahty ( 10.85[ ) becomes 



[d + ur-\Qn) -{d + uT-\Qm) = (10.87) 

which is obviously true since 

[d + ur-\Qn) = {d + ur-\d + <(1) = {d + n)"+"^-i(l). ■ 

Thus, we have found the general l'^*-order extensions 

ut = d{Qn), ift = {aQn-1,^ + Qn-ld){ip). (10.88) 



The remaining case of ( 10.78b| ), /? = 0, leads after a quick calculation with a = to 
the trivial extension 

uj = p = 0; a = b = c = d = e = 0. (10.89) 

It appears that there exist no nontrivial extensions of the Burgers hierarchy given by 
extension operators of order zero. Nevertheless, we have: 

Proposition 10.90. The extended flows 

Ut = d{Qn), ^t = iQnV, n G N, 7 = const, (10.91) 

all coTnmute. 
Proof. We have to show that 

XniQm) - XrniQn) = 0. (10.92) 

But, by (|10.86a|) , 

XniQm) = {d + ur - Qm) {Qn) = Qn+m " QnQm, (10.93) 

and this is symmetric in n, m. I 
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Remark 10.94. The system ( |10.91 ) is a restriction 
of a larger system 



(10.95) 



(10.96) 



from C{f,V'}- The dual to ( 10.91 ) system has the same form ( 10.91 ), with ^'^^'^ = — 7°'^^. 



Whether the system ( [10.84 ) could be extended from C{u,ip} into C{v,'ip} or not, I was 
not able to determine. 

11 The equation Ut + uUx = and its hierarchy 

In the quasiclassical limit, the Burgers hierarchy 

ut = d{Qn), n€N, (11.1) 

becomes 

(X„) nt = a(tz"), nGN. (11.2) 

These are quasilinear equations. Let's determine all their quasilinear extensions, of the 
form 

(X-*) ut = d{n^), ^t = {an{n''-')^ + PnU^-'d){y,). (11.3) 

The linearization of (11.2) yields 

^t = d {nu^-^^) = {n (n"-i)^ + nu''-^d)i^), (11.4) 

so that 

an = n, Pn = n, (11-5) 

while its dual provides 

^t = (nu'^-^d) (if), (11.6) 

so that 

On = 0, (3n= n. (11-7) 

The 9^-terms of the equation ( |1.13|) , 

yield: 



Pmim - l)n - Pnin- l)m = (m - n)PnPm- 



(11.9a) 
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Regrouping, we find 

Pmim - l)(n - (3n) = Pnin - l)(m - Prn)- 
The degenerate case 

/?„ = 0, V n, 
we shall examine later on. The case 

Pn = n, V n, 
corresponds to the linerization ( 11. 5| ) and its dual ( 11.7 ). So, suppose 

Pn + 0, n. 
Rewriting the equation ( |11.9| ) as 

/3m(m-l) /3„(n-l) 



m- i3m 
we see that 



n- I3„ 



(11.9b) 
(11.10) 
(11.11) 
(11.12) 

(11.13) 



(3m{fn — 1) = 9{m — (3^), 9 = const 
9m m 

9 + m—l l+a;(m — 1)' 



oj = const. 



When w = oo, we recover (3jn = ( [LLIOP , and when a; = 0, we obtain (3„ 
Let us now tackle the 5^-terins: 



m 



(11.14) 
dLUD- 



,71—1 t„,m~l 



UnamU^'^ ("'""').. - /3m«nn"^-' (n^-^) 



(11.15) 



I claim that the a„'s are just the div-terms {const /3n}- Indeed, fix n > 1. By Corol- 



lary 9.18 and formula ( 9.21 ), we can make a„ vanish for this particular n. Equation 



(11.15) then becomes 



ar 



(11.16) 



Since /?„ 7^ 0, am must vanish because the RHS of ( 11.16| ) is not a total derivative while 
the LHS is: 



,n — 1 l^.m—X 



K'').. ~ - K"'). (^™"'). = -(^ - 1)("^ - 1)^"+'"-'^' 9^ 0. (11.17) 



Thus, all extensions of ( |ll.2| ) with /3„ 7^ are of the form 

1 + uj\n — 1) ^ 
and their dual 

1 + LO\n — \) 



(11.18) 



(11.19) 
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For general lo, these extensions are rational in n and certainly do not come under quasi- 
classical limit from any Burgers extensions of Section 10. 

It remains to tackle the case /3„ = ( |11.10 ). The equation ( 11.^ ), in the form (11.15), 
becomes, for n,m> 1: 

a^Xn (n"-!) = Q„X„ (n"-i) , (11.20) 

which is 

am{m - l)u"""^nn"~^u^ = an{n - l)u''~'^mu"'~'^u^, (11-21) 

so that 

am{m — l)n = an{n — l)m. (11.22) 

Since n,m > 1, we find 

am('n' — 1) = Ai"^5 // = const, (11.23) 

and we obtain the extensions 

ut = d (u") , (ft = finu^~'^Ux(p. (11.24) 

In this form, not the form ( 11. 3| ), it is applicable also to the case n = 1: 



ut = Ux, (ft = fiu ^Ux^, (11.25) 

which is a rational extension, not a polynomial one. As such, it is not a quasiclassical limit 
of any of the Burgers extensions of Section 10. 

12 Long wave equations 

The 1 + 1 — d system 

ut = uux + hx, ht = {uh)x (12.1) 

for a pair of functions u{x, t) and h[x, t) is the oldest one considered in this paper. Riemann 
analyzed it in 1860 [25]. More recently it was realized that this system [17, 18, 12] (and 
its various generalizations) is a part of an infinite hierarchy of commuting Hamiltonian 
systems of the form 

:)r(ao)(v-)<->' <-* 

with 

H2 = ^ (12.3) 

being the Hamiltonian of the system ( [12. 1[ ). The next flow, ^^3, 

ut = u^Ux + {2uh)x, ht={hu^ + h^)^, (12.4) 
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has the Hamiltonian 
H3 = — + h\. 



(12.5) 



Let us determine which extensions, if any, of the flows (12.1) and (12.4), still commute; 
we assume that the conditions ( |6.1| ) are in force. So, let the extensions of (|12.1| ) and ( |12.4 ) 
be, respectively: 



(ft = {aux + Pud){(p), 

ift = ii'yuux + 5hx) + {fJ.u^ + vh) d){ip), 

where we have utilized the natural grading 

rk{h) = 2rk{u). 
First, 

^ = =^ fi = u = =^ a = j = 6 = 0. 
So, assume 

/3/0. 



By Corollary |9. 18 and formula ( |9.21 ), we can then make a vanish, and this forces 
as well. Thus, our extensions ([L2.6|) become 



ipt = Puifx, ift = {tJ-u^ + i^h) fx- 

Equation ( |1.13 ) then yields: 

2fj,u{uux + hx) + v{uh)x — [5\u Ux + 2{uh)^ 
= f3u [fiu + uh) , — [fiu + vK) Pux- 

Equating the like-terms, we obtain: 

2/i - /? = /i/3, 
1/ - 2/3 = -I//?, 
2fi + u-2p = pu, 



(12.6a) 
(12.6b) 

(12.7) 

(12.8) 

(12.9) 
7 = (5 = 

(12.10) 



or 



^(2 -/?) = /?, 
Kl + /?) = 2/3, 

^ - iy/3 = 0. 

Eliminating fi, v, we get: 



(12.11) 



(12.12a) 
(12.12b) 
(12.12c) 



(12.13a) 
(12.13b) 
(12.13c) 



P 



2/?: 



P 



P 



2-/3 -"1 + /? (2-/3)(l + /3) 
= (2/^2 _ 3/3 + 1) = 2(/3 - 1)(/? - 1/2). 



{1 + /3- 2/3(2-/3)} 



(12.14) 
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Thus, we got two different extensions: 
/3 = 1, fi = i^ = 1 : 

ft = Uifx, 

ipt = {u^ + h)ipx] 

13 = 1/2, /i = 1/3, 1/ = 2/3 : 
ft = -^fx, (pt = -[u + 2h) ifx. 



(12.15a) 
(12.15b) 



(12.16) 



The question now is if these two extensions apply to the whole hierarchy of commuting 
flows — or not. 

We start with the /3 = 1-extension ( |12.15| ) first. Noticing that 



ldH2 2 , , 15^3 

U = — — :^ , U + h 



h du ^ ' h du 

We arrive at the following explanation of formulae ( [12.15| ) . 
Proposition 12.18. The following matrix is Hamiltonian for arbitrary constant a: 

u h if 

1 
d 

(9 





(12.17) 



u 
B= h 



( 



--fx + a-dip 
n n 



1 ^1 

(f \ -rfx + afo- 



(12.19) 



/ 



Proof. Proof is a tedious but straightforward verification of the Hamiltonian criterion 

[14, p. 47] 



B6 [Y^B{X)] = D[B{Y)]B{X) - D[B{X)]B{Y), V X, Y. 



(12.20) 



The matrix B ( 12.19| ) provides a Hamiltonian solution to our extension problem. As 
a bonus, this same matrix automatically supplies an extension of the Dispersive Water 
Waves (DWW) hierarchy [12] 



Ut = UUx + hx+ 9Uxx, 

ht = {uh)x — Ohxx, G = const, 



(12.21a) 
(12.21b) 



because the system ( 12.21| ) is also of the Hamiltonian type ( |12.2| ), with the Hamiltonian 
H2 = ^ + euxh. (12.22) 



Let us now analyse the case (3 = 1/2 ( [L2.16| ). Recall that each of the systems we are 
looking at, the hierarchy of commuting systems including ( |12.1 ), is not just a Hamiltonian 
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system but a in-Hamiltonian system [12]. In particular, the 2 Hamiltonian structure is 
given by the Hamiltonian matrix 



B 



II 



2d 



du 



^ud hd + dh^ 
The corresponding Hamiltonians of the flows (12.1|) and (|12.4) are, respectively 



(12.23) 



^2 = T = 2^^' 
_ hu^ + h' _2^ 
^ 3 3 ^' 

In particular, 

6n2 u 6m n2 + 2/i 



6h 



6h 



(12.24a) 
(12.24b) 

(12.25) 



Comparing the latter expressions with formulae ( [L2.16| ), we arrive at the following expla- 
nations of the latter: 

Proposition 12.26. The following matrix is Hamiltonian for arbitrary constant a: 



h 
du 







B 



II 



(2d 

ud hd + dh —Lpx + cadip \ ■ 

(/3 \ 'fx + OLLpd 



(12.27) 



Proof. Apart from the noninterfering u-part of the matrix (12.27), this matrix is identical 
to the Hamiltonian matrix ( tj.23 ). ■ 



Remark 12.28. The matrix ( |12.27] ) provides the /5 = g-extension for the whole hierarchy 
of dispersiveless long wave systems. What about the full dispersive version, the DWW 
hierarchy, of which ( |12.21| ) is just one flow? Nothing to it. The matrix B ( 12.23| ) is the 
quasiclassical limit of the full 2"'^ DWW Hamiltonian matrix [12] 



B 



-II 



2d du + 2(952 \ 

ud - 2ed^ hd + dh ) 



(12.29) 



which differs from the Hamiltonian matrix B^^ (12.23) just by an extra 2-cocycle 29d^ at 

the u-h place. Therefore, the (/9-extension of i? , as that of B , goes through in exactly 
the same way: 



B 



II 



2d du + 26*52 

ud - 29d'^ hd + dh 
fx + otifd 








(12.30) 
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We now return to our original system ( 12.1 ). Because it's quasilinear, the weight rk{d) 
is undetermined. This forces us to the extension form ( |12.6| ). However, the fuh dispersive 
system ( 12.21| ), of which ( |12.1| ) is the quasiclassical hmit, fixes the weight of d as 

rk{d)=rk{u). (12.31) 

This allows for a far greater apriori form of extensions than ( |12.6| ). Let us examine one 
such extension, of zeroth order. 
We assume it be of the form 

^t = Hn{u, h)ip, (12.32) 

where Hn{u, h) is a polynomical in u and /i, homogeneous of weight 

rk{Hn) = rk{Hn) - 2rk{u) = nrk{u). (12.33) 

Thus, the flows ( |12.1 ) and (12.4) are extended as 

ipt = {aih + a2U^) <f, (fit = {ashu + 04^^) (p, (12.34) 

respectively. The commutativity of the extended flows leads to 

H2=a(h + ^j, (12.35a) 

^3 = a [ 2/in + ^ j , a = const. (12.35b) 

If we notice that formulae ( |12.35| ) can be rewritten as 

we can extend formulae ( |12.35| ) to the whole long wave hierarchy: 

Proposition 12.37. The following matrix is Hamiltonian for arbitrary constant a: 





u 


h 


V 


u 


fo 


d 





h 


d 





—aif 


"P 


Vo 


aip 






B= h\ d -ap ■ (12.38) 



Proof. The h-ip part of -B is a skewsymmetric non-differential 2 by 2 matrix, and so 
is Hamiltonian (and Lie-algebraic). The remaining u-h part represents a 2-cocycle on a 
trivial Lie algebra. ■ 

Again, the Hamiltonian matrix B ( 12.38| ) extends the full dispersive DWW hierarchy, 
not just its quasiclassical limit. 
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Remark 12.39. The Hamiltonian extensions (|12.19| ) and ( 12.38| ) of the Haniiltonian ma- 
trix ( 12.21) , are far from being the only ones possible. For example, if we for a moment 
forget about the homogeneuity requirement and notice that u and h enter into the Hamilto- 
nian matrix ( 12. 2| ) on equal footing, we can interchange u and h in each of the Hamiltonian 
matrices ( |12.ig| ) and (12.38). In the case of the latter, we get the Hamiltonian matrix 



B 



■N 





u 


h 


93 


u 


( 


d 


—a^p 


h 


d 








^ 


\aip 









(12.40) 



However, this Hamiltonian matrix provides an extension that is no longer homogeneous 
for the long wave system ( 12. 1| ) and every other system in the hierarchy, because the 
(/9t-equation generated by the Hamiltonian matrix B^ ( |l 2.401) yields 



rk{dt) = rk{H) — rk{u), 
while the /li-equation provides 

rk{dt) = rk{H) — rk{u) + rk{d) — rk{h), 



(12.41a) 



(12.41b) 



and these weights are not equal. However, the remedy is clear: if is not a good variable, 
but 



If = log{ip) 
is. In the variables {u,h,Lp}, the Hamiltonian matrix B-^ ( 12.40 ) becomes 



(12.42) 



H 




B 



Now the c^j-equation yields 

rk{dt) = rk{H) — rk[u) — rk{<f), 
so that setting 

rk{(p) = rk{h) — rk{d), 



(12.43) 



(12.44) 



(12.45) 



the homogenuity is restored — in the c/5-language. Obviously, formula (12.45) can not be 
transferred back to the original (^-language, because 

ip = exp{ip). (12.46) 

We shall increasingly deal with similar quasi-homogeneous extensions in what follows. 
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Remark 12.47. The hierarchies of long wave equations, both dispersive and dispersive- 
less, allow the reduction 



{h = 0}. 



(12.48) 



Let's see what happens with the dark extensions under this reduction in the zero-dispersion 
case. Here 



dHn 

du 



hnvJ"-^ + O ih^] 



dHn 
dh 



u'^ + Oih), 



and therefore the n flow is 
'dHn 



ut = d 



ht = d 



dh 



= 5(n")+0(/i). 



On the constraint {h = 0} ( 12.48| ) thus becomes 



(12.49) 
(12.50) 

(12.51a) 
(12.51b) 

(12.52) 



which is the flow (|11.2| ). The y?-equation from the Hamiltonian matrix ( 12.191 ) (discarding 
a as inessential) is 



so that on the constraint {h = 0} this becomes 



n — 1 

(ft = nu ¥?j 



(12.53) 



(12.54) 



This is the case Pn = n ( |ll.6D , ( |ll.7D , ( |ll.ll| ). The 2"*^ dark extension, provided by the 
Hamiltonian matrix (12.27), yields: 



u^ = 2d{^^]+du^^^-' 



du 



dh 



h, = ud { ^^] + {hd + dh) ^ ^""^-^ 



(Pt = 'Px 



du 
dh 



dh 



d (nn"-i) + Oih), 
-- Oih), 



<^,n"-^ + 0(/i). 



On the constraint {h = 0} this becomes 

Ut = d{u"-), (/9j = n"~Va;- 

Here /J^ = 1, which corresponds to the case w = 1 of formula ( |11.14 ). 



(12.55a) 
(12.55b) 
(12.55c) 

(12.56) 
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Remark 12.57. The nonlinear Hamiltonian matrix B ( |12.19| ) looks rather unpretensious 
but is in fact very remarkable and stands at the intersection of many interesting fluid 
dynamical and Hamiltonian theories. As a simple illustration, pick an arbitrary function 
E{u), and exchange the variable ip into the variable 



ip = ip + E{u). 



(12.58) 



For Q = 0, the Hamiltonian matrix B (|12.19) becomes, in the new variables {u,h,il)): 



B 



u 
h 




d 



A 



h 

d 

End 



{^ - E)a 

dEu 





where 

A = 1. 
If we now set (without any visible justification) 

A = 0, 
the matrix ( 12.59| ) becomes 



(12.59) 



(12.60) 



(12.61) 



u /O d 
B = hid dEu 

ip \0 Eud 



(12.62) 

This matrix is still Hamiltonian! It can be interpreted as follows. Start with the {u, h)- 
system (|12.2D : 

Now introduce the new variable 

il^ = E{u). 



«'f) 



(12.63) 



(12.64) 



Then 



V'i = EuUt = Eud["-i^]. (12.65) 

The third row of the Hamiltonian matrix B ( |12.62D produces this exact motion equation 
(12.65). However, the third row of the Hamiltonian matrix (12.5S|) yields the same result 



( 12.65 ), since ijj is in fact E{u). Thus, Hamiltonian formalism doesn't discriminate between 
the matrices (12.59) and ( |12.62[ ). (However, the fluid-dynamical considerations show that 
the more complex Hamiltonian matrix ( 12.59 ) is the correct one.) 
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13 The Benney hierarchy 

In 1973 Benney [2] derived the following remarkable 2 + 1 — d generalization of the classical 
1 + 1 — d long wave system ( |12.1 ): 



ut = uux + hx 



ry / rh 

u„ I Uxdy, ht= { udy 
'0 \Jo J X 



(13.1) 



Here h = h{x, t) denotes the height of the fluid free surface, counted from the bottom 
{y = 0}, and u = u{x,y,t) is the horizontal component of the fluid velocity; as always, 
the sign of time t is changed for aesthetic reasons. When u is independent of y, 



0, 



fl3.2) 



the Benney system ( |13.1 ) turns into the classical system (12.1). What makes the Benney 
system so remarkable are the following amazing properties of it that Benney had found: 
1) Introduce the moments of u: 



Ai = Ai{x,t) = / u'dy, i G 
Jo 



Then the system ( [13.1 ) implies: 

Ai^t = Ai+i,x + iAi-iAo,x, i& 



(13.3) 



(13.4) 



2) For each n € Z_|_, the moments system (13.4) has an infinite number of conserved 
densities 



Ha = An 



H, 



Ai, H2 = A2 + Al ... , Hn(^An + Cl[Ao,... ,An-2]. (13.5) 



Later on [17, 18], Manin and myself found that: 

3) (|13.4[) is a Hamiltonian system, in the Hamiltonian structure 



and with the Hamiltonian 

\2\ 



H2/2 = {A2 + Al) /2; 



(13.6) 



(13.7) 



4) All the Hamiltonians Hn ( |13.5 ) commute with respect to the Hamiltonian structure 

(H); 

5) Under the reduction the map ( 13. 2| ): 



Ai = hu^, i G Z+, 



(13.^ 



the Hamiltonian structure ( 13.6 ) collapses into the Hamiltonian structure ( 12.2 ); 

6) For any Hamiltonian H = H{A}, the Hamiltonian system in the moments space 
produced by the Hamiltonian structure (|13.(]| ): 



A. 



n,t 



Y^inA 



n+m-l^ + Sm^n+m-l)(-f^|m); H^^ 



m>0 



5H_ 



(13.9) 
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is lifted up by the moments map ( |13.3 ) from the fohowing system in the physical {u; h}- 
space: 



ut 



E 



;"'"-W|m) 



— u„ 



ry 

Jo 



dy (u 



m—l 



mH\rn), 



^(ylm-l"^^|m)^• 



(13.10a) 
(13.10b) 



We aim in this Section to generalize into the moments language those results of the 
preceding Section which deal with the dispersiveless long wave equations; the full dispersive 
case will be discussed in Section 14. 

We start with the {/5 = l}-case, covered by the Hamiltonian matrix B ( 12.19 ). We'd 
like to extend the Kupershmidt-Manin Hamiltonian matrix B ( p.3.6 ) into a matrix B in 
such a way that the reduction {Ai = hu^} of B reproduces the matrix B ( |12.19| ). How to 
approach this problem? 

The extension motion equation for cp, generated by the matrix B ( |12.19| ), is 



y^t 



-ipa; + aipd- 
n n 



6H_ 

Su 



(13.11) 



Now, under the reduction {Ai = hu"^} ( |13.8| ), 



6u 
5H_ 

~5h 



E 



5Am du 
6H dAm 
6Am dh 



|m! 



E'^i. 



tX^-"' 



\m.i 



(13.12a) 
(13.12b) 



and it is plausible that the 99j-equation ( |13.11| ) comes out of reduction of the equation 

^t = ( ir'fi^ + "V'C^^r ) yZ mA.m-iHi^. (13.13) 

V^o ^0/ — 

I say "plausible" because the RHS's of formulae (13.12) are hugely nonunique, so the 
equation ( 13.13| ) should be considered at the moment as a guess in need of verification. 

Proposition 13.14. The following matrix B is Hamiltonian for arbitrary constant a: 



Br, 



nAn+m-id + dmAn+m-1, n,me 



B^m = ifx + aipd) 



mA, 



m— 1 



A. 



B 



uAn- 



An 



-{adip- ipx), 



B^^ = 0. 



(13.15a) 
(13.15b) 
(13.15c) 



Proof. We need to verify the identit y (|12.2C| ). A part of this verification, corresponding 
to the submatrix {Bnm} ( 13. 6| ) in B ( 13.15]) , can be omitted, since {Bnm} is known to be 
Hamiltonian. The remaining verification is long, mind-numbing, not illuminating, and so 
is better left out. I had done it, and the identity ( 12.20| ) is true. ■ 
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Example 13.16. For the Benney system ( |13.7 ), 
1 



Hi, 



so we get 



-,54 + AoC, 



^t = i'fx + aifd) 



An 



[ip^ + aipd) i- udyj. 



(13.17) 

(13.18a) 
(13.18b) 



This equation becomes ipt = mpx ([12. 15a) when 



Next, let's take up the Hamiltonian matrix B ( |12.38 ), which, in the variable 

^ = log(<^), (13.19) 

takes the form 



(13.20) 




Since the c/?t-equation 
5H 



Lpt = a- 



Sh 



can be rewritten, by formula ( 13.12b| ), as 



we arrive at the following matrix B: 



a 



nAr, 



nm — "-^n+m+ld + OmAn+rn-l, n, m, € 



A A 

^0 ^0 



e. 



Vf 



0. 



(13.21) 



(13.22) 



(13.23a) 
(13.23b) 
(13.23c) 



Proposition 13.24. The matrix B ( 13.25 ) is not Hamiltonian for a y^ 0. 



Proof. The (p-entiy of the criterion (12.20), vanishes on the LHS and doesn't on the 
RHS. ■ 
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Thus, the extension ( 12.35a| ) 



(13.25) 



and similar ones for the higher flows, remain unhfted into the moments space. The source 
of the difficulty is easy to pinpoint: it is the RHS of the formula ( 13.12b| ), where n™ was 
represented as Am/Ao (when Uy = 0), but could have been equally well represented as 
Am+i/Ai or just about any other rational function of {A}. 

We would encounter exactly the same problem with the Hamiltonian matrix 13 (12.27), 
since there 



^t = ifx + aipd) {5H5h) 



(13.26) 



The (/?- independent part of the matrix B ( tl2.27| ) , the Hamiltonian matrix B ( |12.23D , 
comes out of the second Hamiltonian structure of the Benny hierarchy [10, p. 373]. It is 
equation ( 13.26| ), like the equation ( 13.21 ) before it, that resists momentous understanding. 



The same problem appears in the fully dispersive case, exemplifled by the Hamiltonian 
-II 



matrix B ( p^3o|) . 



The last remaining Hamiltonian matrix of Section 12, B^ ( 12.43 ), has no non-zero 
5H / 6h-e,iAx\es in the (p-iaw; there 



'^t 



a- 



5H_ 

5u ' 



(13.27) 



and by analogy with formula ( |13.13| ) , we can utilize formula ( |13.12a ) to guess the following 
generalization of the matrix B^ ( |12.43| ): 



■H 
nm 

•H 

ipm 

'H_ . 



: nAn+m+d + dmAn+m-1, n,m e 
-- amAm-i, BH^ = -anAn-i, 
0. 



(13.28a) 
(13.28b) 
(13.28c) 



Proposition 13.29. The matrix B ( 13.28 ) is Hamiltonian for any constant a. 



We defer a Proof until Section 15, as the matrix B ( 13.28| ) hides some interesting 
mathematics behind it. At the moment let's look at how this matrix could have been 
derived (and in fact was) by considerations entirely within the moments space. 

The Kupershmidt-Manin Hamiltonian matrix B ( 13. 6| ) has two distinguished Hamil- 
tonians: V = Ai, a momentum: 



ax 
and H = aAo, a Casimir: 

^aAo = 0. 

If we extend the momentum flow X-p by the Casimir oAq : 

(ft = aAo 



ieZ+, 
(<^ ft = aAo(p), 



(13.30) 

(13.31) 

(13.32a) 
(13.32b) 
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and ask for such an extension of the arbitrary flow ( |13.9| ) with an x-independent Haniil- 
tonian H, of the form 



that it commutes with the flow ( 13.32| ), we shah get 



\ "1 / X 



so that 

7{H) = y^aniAm-iHijn, 



(13.33) 



(13.34) 



(13.35) 



and we thus recover the ;B^^-entry ( 13.28b| ) of the matrix B^ ( |13.28 ). 



14 The KP hierarchy 

The Benney hierarchy of Section 13 is the quasiclassical limit of the KP hierarchy [27, 4] 



Lt = [P+,L], 

oo 
i=0 

The (1**) Hamiltonian structure Bnm (|13.6| ) of the Benney hierarchy, 



Bnm = nAn+m-id + dmAn+m-i-, n,me 



is the quasiclassical limit of the (1** ) Hamiltonian structure of the KP hierarchy 



B„ 



E 



n 



/^ 



-'^n+m—fiC/ 



-dy 



m 



Ai 



A. 



n+m—fj, 



n,m G Z+. 



(14.1) 
(14.2) 
(14.3) 

(14.4) 
(14.5) 



We aim in this Section to find dispersive analogs of the two zero-disperson results of 
the preceding Section, the Hamiltonian extensions B ( 13.15|) and B^ ( 13.28|) . 

We begin with the Hamiltonian matrix B (13.15), an 1-dimensional extension of the 



Hamiltonian matrix B ( 14.4 ). Recall that the latter Hamiltonian matrix is connected via 
the reduction map 



Ai = hu\ i G 



(14.6) 



with the Hamiltonian matrix b ( 12.2 ) 
u h 

h\d / 



(14.7) 
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Similarly [12, p. 66], the full KP Hamiltonian matrix B ( |14.5| ) is connected to this same 
Hamiltonian matrix b ( |14.7 ) via the differential reduction map 



Ai = hQi{u), ieZ+, Qi{u) = {d + uy{l), 



(14i 



a dispersive analog of the zero-dispersion reduction map ( 14.6| ). Now, the motion equation 
for (f generated by the Hamiltonian matrix {b ( |14.7 ) extended into the Hamiltonian matrix 
(|12J[9|) }, is: 



(ft = {(px + a^ojj^-r^, 
and under the reduction map ( |14.8| ), 



6H 
6u 



E 



DArn 

Du 



6H 



DQr, 



By formula (|10.28a| ). 



DQ„ 



Du 



E 



m 



^m—s(^ 



A'=Y.yr-i":) 



^ s>0 ^ ^ 



DQr 

Du 

s>0 



^-E(-^)' 

s>0 



m 

s + 1 



m 

s + 



s + 1 

■^va—l- 



-, 1 ^m-l-s(-ff|r 



We thus arrive at the equation 



Hence, a plausible extension of the KP matrix B ( 14. 5| ) is: 






n+m—fj, 



Wm—s 



n,m G 



m— 1— S; 



Bnip = ^[ cj + 1 j ^n-l-sO'' — {adip-ipx), B^^ = 0. 



(14.9) 



(14.10) 



(14.11a) 



(14.11b) 



(14.12) 



(14.13) 



(14.14a) 
(14.14b) 
(14.14c) 



Dark Equations 



421 



This matrix B has the following 2 properties: 

1) Its quasiclassical limit is the Hamiltonian matrix B ( |13.15| ); 

2) Under the reduction Ai = hQi{u) ( 14.8 ), the matrix B ( |14.14| ) reduces to the Hamil- 
tonian matrix ( |12.19| ). 



Proposition 14.15. The matrix B (14-14) ^^ Hamiltonian for any constant a. 



Proof. We have to verify the identity (12.20) for the matrix B: 
B5[Y^B{X)] = D[B{Y)]B{X) - D[B{X)]B{Y), 



(14.16) 



with X = {X,u) and Y = iX,v). 

To say that this verification is long, mind-numbing, and not illuminating, would be to 
indulge in coy understatements. It is all that, but in the end one arrives at some non- 
evident identities to be checked, so I provide below a guide to the verification sequence. 

1) Introduce the variable (p = log(c/9) instead of ip. The operators {p^x + otipd) in B^pm 
and {adif — tpx) in Bn,fi become {ad + (px) and {ad — <Px), respectively; 

2) The variables u and v, after grouping and cancellations, should enter only in combi- 
nations 



U 



aui 



ipxU 



An 



V 



avx 



VxV 



An 



(14.17) 



3) In checking the c^-entry of the vectors in ( 14.16| ), the U- and V-terms cancel out, and 
the remaining expressions can be divided out from the left by the operator {ad + ipx)-^- 
What remains is a trilinear differential identity m. X = Xm, Y = Yn, A = Al, for all fixed 
{m,n,L). To prove it, pass to the symbols by substituting 



d « d/dz, X^e" 



Y 



oV^ 



A 



(14.18) 



and then use Newton's binomial summation repeatedly; 

4) In the ^Tv-entry of the vectors in ( 14.16| ), the terms bilinear in {u,v) check out; 
the terms bilinear in {X, Y) are accounted for by the known Hamiltonian character of the 
submatrix B (|14.l 



of B (14.14); the remaining terms are bilinear in {X,V) and (Y,U), 
and due to the skewsymmetry of B, only the (X, V)-entries need to be checked out. After 
repeated cancellations and simplifications, one arrives at the identity 



u-\-k=L 


'( N\ 
\ ^ J 


Ad"- 


-(" 


I ^^{-QfA 


1 

X 

\ 








= - E 

u+k=L 


V 


) 


) - f "" ^ {-dyx 
\ ^ ) 


1 

A 
\ 


' N + m-i^ \ (fc) 
, k+1 J \ 


+ ' 


,,— T ^ 


N ^ 

c+lj 


v(^-) 




^) _ 


^"){-dr{^ 


IX) 


5 



, (14.19) 



for all fixed N, m, L. Passing to symbols, multiplying by vt ~^^, and summing on L, we 
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arrive at the identity 

-[1 + r(x + v)f [(1 + Tvr - 1] + [1 - r(a + X + v)Y 



1 + 



TV 



(1+Tf)^+™ 



1 — T{a + X + V) 

TX 



1} = {1+TX 



,N 



1 + 



1 + TV 



N 



+ {1 + TV 



.N+m 



[l-T{a + x + v)Y 
Tia + X + v) 



1 + TV 



[1 - T{a + x)r [(1 + Tv)^ - 1] + (1 + 



TX 



,N 



1 + 



TV 



1 + TX 



N 



1 



(14.20) 



which is obviously true. 



Our next extended Hamiltonian matrix is B ( |13.2§| ). To find its dispersive analog we 
use the motion equation ( 13.27 ) 



(ft = «■ 



6H_ 

6u 



and formula ( |14.21| ). We find 



ipt = a 



gH 



Ei-^y 



m 
s + l 

The resulting matrix B is 

2H 



A. 



m—l—s 



{H\m)- 



n 



fJ' 



■^n+m—^CI 



-dy 



m 



fJ' 



A, 



n+m—^ 



B- = a 



E(-^)^ 



m 

s + l 



im— 1— S) 



B 



•H 



-a 



n 
s + l 



A,_l_,a^ B, 



•'H_ 



0. 



(14.21a) 
(14.21b) 

(14.22a) 
(14.22b) 
(14.22c) 



This matrix B (14.22), like the matrix B ( 14.14|) bef ore it, has the following 2 properties: 

1) Its quasiclassical limit is the matrix B^ ( [L3.28| ); 

2) Under the reduction Ai = hQi[u) ( |14.8D , the matrix B reduces to the Hamiltonian 
matrix B'^ ( |12.43| ). 

Proposition 14.23. The matrix 8^ ( 14-2^ ) is Hamiltonian for any constant a. 



Since this matrix is linear in the variables {A}, there is a Lie algebra lurking behind it. 
We examine the mathematics responsible for this matrix to be Hamiltonian in the next 
Section. 

At the moment, let us derive the matrix B^ ( 14.22 ) in purely {A}-language, like we 
did at the end of Section 13 with the quasiclassical limit B (|13.28D of this matrix. 

First, H = Aq is a Casimir of the Hamiltonian matrix B ( 14.5 ), because 



b™=e(^ )■*-.*'= e(, TO-'"-'-'**''" 



M>0 

Xao (An 



0, ViVG 



(14.24) 
(14.25) 
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Next, V = Ai is a momentum: 

^TVi = ^ ( ] AN+i-i.df' - An+1 + OAn 



E 



N 

/U+1 



AN-ad^+^ + dAN 



(14.26) 



=> XaMn) = d{AN), V iV G Z+. (14.27a) 

Therefore, if we extend the flow X-p by the equation 

(ft = aAo, (14.27b) 

and every other flow Xh with x-independent Hamiltonian H, by the equation 

'ft =nH), (14.28) 

and require that the thus extended flows still commute with the flow ( |14.27| ), we find: 

XaAUH)] = [7{H)], = XniaAo) = a J]5o™(/7|^) 

[by adjoint of (|l4.24|)] v^ / r^^„4_^ , / m 



-a 



^(-a)^+i^„_i_ 



m,ii 



Ai + 1 



i^i 



\mj 






m— 1— /x \J^\m) 



(14.29) 
(14.30) 



Substituting this in ( 14.28 ), we recover precisely the equation ( 14.21b|) and hence the 
matrix 5^ ( pX2^ ). 



Remark 14.31. The KP Hamiltonian matrix B ( |14.5| ), under the differential reduction 
Ai = hQi{u) ( 14.^ ), reduces to the 1 + 1 — d Hamiltonian matrix h (14.7); the KP hierarchy 
is thereby reduced to the dispersive DWW hierarchy of Section 12. Similar conclusion 
applies to the quasiclassical limit of this picture, with the Benney hierarchy being reduced 
to the zero-dispersion long wave hierarchy of Section 12. However, the Benney hierarchy 
also comes out of a genuine 2 + 1 — d hierarchy under the moments map Ai = L u^dy 
( 13. 3| ). There must exist a common picture uniting both these properties, dispersion and 
2 + 1 — d, something like the map 



Ai 



Qi{u)dy, 



(14.32) 



reproducing the KP hierarchy out of some 2 + 1 — d dispersive hierarchy. None has been 
found so far, and this remains one of the most important unsolved problems in the theory 
of integrable systems, the true Holy Grail. 
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Remark 14.33. In this Section we have constructed two 1-component Haniiltonian ex- 
tensions of the KP hierarchy, those given by the Haniiltonian matrices B (14.14) and B^ 



(14.22). These extensions are different from still another l-component Haniiltonian ex- 



tension of the KP hierarchy, the one provided by the mKP hierarchy ([15] and [16, § 6.2]). 
Formula (6.2.3) in [16] shows that that third extension is in fact of the extended Lax 
type of Section 5; the unexpected fact proven in [16] is that that Lax-type extension is 
Hamiltonian — we have found nothing similar in this paper, even though the Hamiltonian 
property of extended Lax systems could not be ruled out with certainty. 

15 One-dimensional linear extensions 
of linear Hamiltonian matrices 



In Sections 13 and 14 we derived two conjecturally Hamiltonian matrices, B (|13.28 ) and 
B'^ ( [14:221) : 



B'^^ = nAn+rn~id + dmAn+m-i, n,meZ+, (15.1a) 

B'^^ = amAm-i, B1^^ = -anAn-i, (15.1b) 

B'^^ = 0, (15.1c) 






I -^n+m—fj.'J V C*) I I An-\-rn—fi 



^15. 2a) 



B^m = « E(-^)' i .Z^) ^—1- (1^-2^) 



^^^ = -«E ,l^ ^4n-l-.5^ B'^^ = 0. (15.2c) 



The 1** of these matrices, B^ ( [L5.1 ), is supposed to be the quasiclassical limit of the 2 



nd 



B^ ( 15. 2| ), although some work is required to make this supposition precise. 



Each of these two matrices has the same general structure: a linear Hamiltonian ma- 
trix; extended linearly by one extra component; with extended entries involving only old 
variables; and with no self-interaction (i.e., B^^ = 0.) 

Since a linear Hamiltonian matrix B = B{q} corresponds uniquely to (the dual space 
of) a Lie algebra, say Q, via the relation (see [14]) 

X^B{Y) r.q^[X,Y], VX,y, (15.3) 

where ^ stands for equality modulo trivial elements ("divergencies"), let us see exactly 
what kind of Lie algebras we get with these types of one-dimensional extensions. 

So, suppose we start with a Lie algebra Q, which we assume for the time being to be 
finite-dimensional-like: this means that no functional operations — like derivations and 
shifts — are involved in the commutators — both in Q and in the 1-dimensional extension 
G of Q. We choose a basis (ej,e) in Q, with (e^) being a basis in Q. Denote by q the 
coordinates on G* , so that the Hamiltonian matrix b = b{G) attached to Q* has the form 
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where < c^- > are the structure constants of Q in the basis (cj): 



ij k 

Thus, our extended matrix h = b{Q) looks hke 



qk 



Y,X^e,.Y.^'e, 



(15.5) 



9i 



b 



f 



EOfQk 

k 



Eokk 







(15.6) 



where 0^ are some constants. Now, by formula ( |15.3| ), the commutator in Q associated 
with the matrix b ( |15.6D can be extracted from the following calculation 

j 

= E 9M E 4^'^* + E ^i K' - ^^') ^ • (15.7) 



V 



Therefore, if we denote hy 9 : Q ^ G the operator acting on Q by the rule 



(15.^ 



the commutator in G has the form 



X 

u 



Y 

V 



[X,Y] + u0{Y)-ve{X) 




The Jacobi identity for the commutator ( 15. 9| ) is equivalent to the equality 

[e{x), Y] - [eiY),x] = e{[x, y]), ^x^ygg. 



(15.9) 



(15.10) 



Thus, the matrix b ( |15.6| ) is Hamiltonian iff the operator 9 : G ^ G is a derivation of G- 
In other words, there is a one-to-one correspondence between derivations of a Lie algebra 
G and linear one-dimensional extensions of the Hamiltonian matrix b{G) attached to G- 

Although it may be not immediately obvious, this natural one-to-one correspondence 
between derivations and one-dimensional extensions is a purely finite-dimensional-like phe- 
nomenon which doesn't carry over to functional Lie algebras. We shall see about this 



breakdown presently. Notice that both Lie algebras G entering formulae ( 15.1a ) and 



( 15.2a ) are differential, not finite-dimensional-like, so our analysis above doesn't apply to 
the two matrices we are interested in, B^ ( p.5.1 ) and B^ ( 15. 2| ). Thus, we need to analyze 
each one of these matrices directly. 
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For the matrix B , (|l5.lj ) formula ( |15.3 ) yields: 

[{X,u),{Y,v)f = Y, (riX"y™ - my"^X^) (15.11a) 

n+m=k+l 

+ aik + I) (uY^+^ - vX^+^y (15.11b) 

To get a better handle on this commutator, let's pass on to the generating functions, 
considering instead of the vector X = (X"(x)) the function X = X{x,p) = Yl X'^p'^. 

n=0 

Formula ( 15.11 ) then becomes 

oo 

\x,u),(Y,v)]=Y,p'[{X,u),{Y,vt 

fc=0 

= Y, [nXy^-'^ (y>™)^ - mY"'p'^ (^ V) J 
nm 

+ a ^ [n(A; + l)y'=+ V - v{k + 1)X'=+ V" 

k 

= [Xp% - X,Yp^ + a (uYp - vXp^ (15.12) 

= [x,y} + a (Ix,- j v\ + I- j u,y\'\ 

= Ix -a u,Y -a v\ , (15.13) 

where 

{a, b} = apbx - a^bp (15.14) 

is the ususal Poisson bracket on the (3;,p)-plane. Thus, finally, our commutator is 

7 J . (15.15) 



X 

u 



Y 

V 



X — a J u,Y — a J ' 




The Jacobi identity for this commutator follows at once from that for the Poisson bracket 
(15.14) on the (x,p)-plane. 

Formula (|15.15| ) suggests the following general 

Proposition 15.16. Let G be a Lie algebra over (functional) ring TZ and let 7i d Q be 
an abelian subalgebra of Q. Let O : TZ ^ 7i be a linear operator. Then the following 
commutator defines a Lie algebra structure on Q + TZ: 



X 

u 



Y 

V 



[X + 0{u),Y + 0{v) 




(15.17) 



Proof is obvious. 

Formula ( 15.15| ) is the particular case 



[0{u)f 



aJg ' "" 



(15.18) 
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of the Proposition 15.16| , for formula ( 15.11a| ) shows that 7i = 7^{(5o} is an abehan 
subalgebra of Q (|15.11a| ). (More generally, if M is a manifold and Q = C'^{T*M), the 
Lie algebra of smooth functions with respect to the standard Poisson bracket on T*M, we 
can take n = C°°{M).) 

The 2"'^ case, B^ ( |15.2 ), is rather similar. First, for the Lie algebra Q itself, formula 
(|15.2a|) yields: 



n+m=k+fi 






(15.19) 



Passing to the generating functions X = Y^X^p'^, etc., we can rewrite formula ( |15.19| ) as 

k 



x,y]=^|xv--(^) 



M / - 1 (ym^m^M _ ym^m-/. ( ru \ ^^n^n^i.^,) . (15 20a) 



A^ 



XoY -Y oX, 



(15.20b) 



where o is the associative multiplication of (symbols of) differential operators: 



xo Y = y- 



d_ 
dp 



X 



d_ 

dx 



Y 



(15.21) 



Second, for the extended algebra Q we get from formulae (|15.2| ) and ( |15.3 ): 



[{X,u),{Y,v)f = [old[X,Yf} 



a 



y^ , k + S + 1 \ /^k+s+l^is) _ Yk+s+l^is)\ 



(15.22a) 
(15.22b) 



Passing to the generating functions and noticing that 



y^ I k + S + l j j^k+s+lpk^{s) 



k\s 



S + l 



X, / V 



Xol V 



v]oX, (15.23) 



we can rewrite formulae (15.22) as 

/ \X,Y 



u \ V 



a 



X,fv +a Y,Ju 




X — a f u,Y — a f V 




(15.24) 



This formula is identical in form to formula ( |15.15| ) , and therefore we get a Lie algebra for 
the same reason as before, with the map O : TZ ^ Ti oi Proposition 15.16| being exactly 
the same ( 15.18 ). 

Thus, both matrices B ( 15.1 ) and B ( 15. 2| ) are Hamiltonian. 
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Let us notice that in a finite-dimensional- like framework, formula ( 15.17 ) reduces to 
the following particular case of formula (|15.9| ) : 

[0{u), Y] = [uO{l),Y] = u[0{l), Y] (15.25a) 

^ e{Y) = [0{1),Y] = ado(i)(y). (15.25b) 

We now address the problem of linear one-dimensional extensions in general. Let Q be 
a Lie algebra with a fixed basis, and let B = B(Q) be the associated Hamiltonian matrix 
(|15.3| ), so that 

Y,Qk[X,Yf ~ X'B{Y) = Y,X^B,, {Y^) . (15.26) 

k ij 

Let the extended matrix B be of the general form 
Qj ^ 
5 = ^M ^'^ ~^i I , (15.27) 

^\nj J 

where n^'s are some (/^-independent operators linearly- dependent upon the g's. Then 

{X, uYB{Y, v) = X^B{Y) +uY^nj (Y^) -vY^ ni{X') (15.28) 

J i 

~ Y, qk[X, Yf + YQk {mHu, Y) - M\v, X)) , (15.29) 

k k 

where AA.^{u,Y) are some bilinear operators satisfying the relations 

Y,(lkM''{n,Y)^uY,^i{Y') (15.30) 

k j 

Thus, the multiplication in an extended Lie algebra Q D Q corresponding to the matrix B 
([1537D is 



X\ Y 



[X,Y]+M{u,Y)-M{v,X) ^ ^^^^^^ 



Set 

MuiY)=M{u,Y). (15.32) 

Then, for each n, A4u is a linear map of G into Q, and the finite-dimensional-like formula 
( 15. 9| ) is the particular case 

Mu{Y)=u9{Y). (15.33) 

The Jabobi identity for the commutator ( |15.31| ) amounts to: 

= {[[X,Y] + M{u,Y)-M{v,X), Z]-M{w,[X,Y] + M{u,Y)-M{u,X))} + c.p. 
= {[MUX),Y] - [MUY),X]-MU[X,Y]} + c.p. (15.34a) 

+ {-MuM^iZ) + M^Mu{Z)} + c.p., (15.34b) 
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so that the matrix B ( 15.27| ) is Haniiltonian iff the commutator ( 15.3l| ) defines a Lie 
algebra iff 



Mui[X,Y]) = [MuiX),Y] + [X,MuiY)], V X,Y,u, 
MuMv = MvMu, y u,v. 



(15.35a) 
(15.35b) 



Thus, for each fixed u, Mu : ^ — *• ^ is a derivation of G, but as u varies, the deriva- 
tions {^Au} sweep a commutative family. (Formula (|15.17 ) is the case Aiu = ado(„), 
0{u) C TC, 7i being abelian.) We see that, in general, linear 1-dimensional extensions of 
linear Hamiltonian matrices correspond to much much more than just derivations of Lie 
algebras; instead, one must have a commuting one-parameter family of such derivations, 
the parameter ranging over the ring(s) over which our Lie algebras are free modules. 

Notice that if the family Mu satisfies the defining relations ( |15.35| ) and a is a constant 
then the family 



Mu = aMu 



(15.36) 



also satisfies the relations ( 15.351 ). This accounts for the appearance of the free parameter 
a in formulae (15.1) and (15. 2|). 



16 Nonlinear derivations of Lie algebras 

In working out various scalar extensions 



iia.i) 



we have met in this paper a few such where the extension operators A^ are of order zero, 
i.e., functions: 



ut = X^, (pt = ipAiX). 
In this case, and in this case only, one can pass to the variable 

<f = log{ip), 
so that the extensions ( |16.2[) become 

ut = X-, ^t = A{X). 
The examples of this sort include the system ( |10.91 ): 

ut = d{Qn), ft = aQnf, n G N, a = const, 
the system (pT2^ : 

Ut = d (u") , ft = a.nu^~ Uxf, n e N, a = const. 



the system (|12.38D : 



5H 

ft = oif^Ti O! = const, 
oh 



(16.2) 
(16.3) 
(16.4) 
(16.5) 
(16.6) 

(16.7a) 
(16.7b) 
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the system (|12.40D : 

ut = aif——, a = const, (16.8b) 

ou 

and the Hamiltonian matrices B^ ( 13.28|) and B^ (|14.22 ) studied in the preceding Section. 

Looking the above examples over, we may notice that none of them is of a universal 
character (C) of Section 1: the first two, (16.5) and (16.6), are of the type (B), and the last 
four are of the Hamiltonian type (D). As the analysis in the preceding Section suggests, 
this absence of ( (7)-type examples is surely a consequence of infinite-dimemsionality — in 
other words, of working with PDE's rather that ODE's, or with Field Theory rather than 
Mechanics. So let us consider the finite-dimensional-ODE-Mechanics case now. 

We start with a few examples, of the following sort: given a finite-dimensional Lie 
algebra Q whose elements are represented as vector fields on a finite-dimensional space V, 
find a linear function A : G ^^ C°°(y) satisfying the equation ( 1.11 ): 

Ai[X,Y]) = XiAiY))-Y{A{X)), yX,Yeg, (16.9) 

where X is the vector field on C°°(y) corresponding to the element X £ Q. 
Our first example is ^ = s0(3), with the basis (cj), i € Z3, satisfying 

[ei,ei+i] = ei+2, i € ^3- (16.10) 

Set 

d 
ei = zdy-ydz, €2 = xd^ - zd^, e3 = ydx-xdy, 5^ = -^,... (16.11) 

ox 

Let O : so{3) -^ so{3) be the following automorphism: 

0{ei) = e,+i, i G Z3. (16.12) 

Extend O to act on C[x, y, z] in such way that 

OiXif)) = OiX)iOif)), yXGG, yfeC[x,y,z]: (16.13) 

0{{x;y;z})={y;z;x}. (16.14) 

If we demand that solutions of the equation (|16.9| ) be O-invariant — and if we don't, we 
get too many of them — we arrive at the single equation 

ei{0{A))-e2i^) = 0-\A), A:=Aiei), (16.15) 

which is 

d d\ f d d\ 

z-Q- -y-Q^j {Mv^z.x)) - ( x— - z-Q^\ {A{x,y,z)) = A{z,x,y). (16.16) 
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A quick inspection shows that 

deg (^) = 1 ^ A = const {y - z), (16.17a) 

deg(^) = 2 => A = const{y - z){x + y + z), (16.17b) 

and that 

A{x,y,z) = {y-z)f{x + y + z) (16.18) 

is a solution of the equation ( |16.16| ) for any function /. It seems hkely that all polynomial 
solutions of the equation ( [16.16| ) are of the form ( 16.18| ) . 



Our next example is the 3-dimensional Lie algebra s£2 acting on C°° (M^) : 

X = dt, Y = tdt, Z = t^dt, (16.19a) 

[X, Y] = X, [Z, Y] = -Z, [X, Z] = 2Y. (16.19b) 

The solution space of the equation ( |16.9| ) is infinite-dimensional, with a basis 

^(X)=r, AiY) = e+\ A{Z)=e+^, Va. (16.20) 

More generally, consider the centerless Virasoro algebra 

[en,em] = (n- m)en+m, n,meZ, (16.21) 

of which the Lie algebra (|l6.19b| ) above is the subalgebra X = ei, Y = eo, Z = e_i. The 
Lie algebra ( |16.21| ) acts on C[t, i~^] as 

en = t^-'^dt, (16.22) 



and a basis of the solution space of the equation (16.9) is 



A{en) = t""", neZ, V a. (16.23) 

Having the multitude of examples at our disposal, we can now interprete them. 
If (/ is a (finite-dimensional) Lie algebra with a basis (ei), 

h,e,] = 2Z4efc, (16.24) 

k 

then there exists a natural representation of Q by vector fields on C°°(t/*): 

e, = Y,ctqk^, (16.25) 

ks ^^^ 

based on the identification of any element X ^ Q with a linear function ( , X) on Q* 
whose Hamiltonian vector field is ( |16.25| ). This is the origin of formula ( |16.11| ) for Q = 
so(3). A linear function on G* is the same as an element of Q; therefore, a linear map 
A : Q ^f {linear functions on Q*} is just a linear operator A : Q ^ G; the equation (16.9) 



for such A simply says that ^ is a derivation of Q. (The solution (|16.17a) corresponds 
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to A{X) = const [cj + 62 + 63, X].) We can interprete the equation ( 16.9| ) for the case 
A: Q ^f C°°{G*) as defining a nonlinear derivation of Q. 

Moreover, the Hamiltonian structure on C°°{G*) is not too important. Let p : G ^ 
End (W) be a representation of ^ on a vector space W. Now, any hnear operator a :U ^ U 
on a vector space U defines a hnear vector field a on C°°{U*), by the rule 

a{{u, )) = {a{u), ), \/ueU, (16.26) 



where (n, ) is a linear function on U* . Thus, formula (16.26) defines the action of the 



vector field a on linear functions, and the Leibniz rule then extends the action of a onto 



the whole C°°{U*). Formula ( 16.25 ) is a particular case {U = Q, a = [e^, ]} of formula 



( 16.26 ). Therefore, in general, if we set 

X = ^), yXGG, (16.27) 



the equation (16. 6|) is an equation for an unknown linear map A : Q ^ C°°{W*). If 



Im(^) C W C C^iW*), the equation (|16.6D becomes 



Ai[X,Y]) = piX)iAiY))-piY){A{X)), yX,Yeg. (16.28) 

Such a linear map A : G ^ W can be thought of as a relative derivation of Q with the values 
in a (/-module. We get the usual derivation when the representation p : ^ — > End (W) is 
just the adjoint representation. (An analog of an interior derivation of G, A{X) = [X,a] 
for a fixed a E ^, is the following relative interior derivation: 

A{X) = p{X).v, fixed veW. (16.29) 



Equality ( 16. 9| ) then follows from the representation property of p.) The general equation 



( 16. 9| ) for the linear map A : Q ^>- C°°{yV*) describes then nonlinear relative derivations. 



Formula ( 16. 9| ), rewritten as 



X{A{Y)) - Y{AiX)) - Ai[X, Y]) = (16.30) 

suggests that one can get this equation as a 2-cocycle condition in a suitably constructed 
complex. This is left as an Exercise to the reader. 

17 Linear extensions as Symbols of nonlinear ones 

In this Section we return to the unfinished business of the mysterious and the strange 
extensions of Section 6. We prove that such extensions exist for the whole KdV and 
mKdV hierarchies. 

We start with the mysterious extension ( |6.32| ): 

ut = d {3u^ + u^^) , ^t = {-d{6u + 2d^)){v). (17.1) 

This can be recognized as (^-linearized part of the Hirota-Satsuma equation [5]. The 
following explanation is inspired by the Wilson succinct analysis [28]. 
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Let C be an even-order scalar Lax operator: 

27V-2 

£ = a^^ + Y, Uid\ (17.2) 



The commuting hierarchy of Lax equations 

Ct = [V+, C], ord(P) = 2n - 1, n G N, (17.3) 

is compatible with the self-adjointness constraint 
C^ = C. 

So let C be self-adjoint: 

AT-l 

c = d^^ + Y^ {u,d^' + d'^'Ui) . (17.4) 

j=0 

Set 

L = a2 + u, u = lc/^_i. (17.5) 

Then L is self- adjoint, and so are L^ and 

Af-2 

C-L''=Y, ifid^' + d^'fi) , some /^'s. (17.6) 

j=0 

Since 

Lt = [P+,L] => (L^)^=[P+,L^], (17.7) 

the motion equations ( |17.3| ): 

Ct = [V+, C], C) = C, V^ = -V, oTd{V) = 2n + 1, (17.8) 

have the form 

ut = X-{u} + 0{f), (17.9a) 

n,t = 0.,{f) + O (f) , i = 0,...,N-2, (17.9b) 

where ut = X~{u} is the n KdV equation, and Cj's are some linear differential operators 
whose coefficients depend upon u. 

We now take the /-Symbol of the system ( |17.9| ) by setting 

f, = e^i, i = 0,...,N-2, (17.10) 

and then letting e — > 0. We obtain: 

Ut = X-{u}, (17.11a) 

^^,t = Oii^), i = 0,...,N-2. (17.11b) 
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This is an (A^ — l)-component extension of the KdV hierarchy, which should be called the 
dual-to-mysterious {N — l)-extension. This name is justified by the following observation: 
Take iV = 2, so that iV - 1 = 1, ord(£) = 4: 

c = d^ + Uod^ + d^Uo + Vi = {d^ + uf + /. (17.12) 

With oTd{V) = 3,V = 4/:3/4. 

r+ = id^ + 3ud + 3du, (17.13) 

we find: 

Ct = LLt + LtL + ft= {d^ + u)ut + ut {d^ + u) + ft = utd^ + d'^ut + 2uut + ft 

= [V+,L^ + f] = [T+,L]L + L[V+,L] + [P+, f] {with ut = [P+, L] = Quu^ + u^^^} 

= utL + Lut + 4.{fxxx + 3fxxd + 3fxd^) + 6ufx 

= utd^ + d\t + 2uut + 6fxd^ + d^6fx - 6(/,,, + 2fxxd) + ifxxx + 12/^.5 + 6ufx 

= {ut + 6fx)d'^ + d^{ut + 6fx) + 2uut + Qufx - 2fxxx 

=^ ut = ut + dfx = Quux + Uxxx + Qfx, (17.14a) 

ft = -2uut + 2uut + Qufx - 2fxxx = -Qufx - 2fxxx- (17.14b) 

Taking the /-Symbol of the above, by setting 

/ = £¥5 

and letting e go, we find 

Ut = d {3u^ + Uxx) , vt = -{6ud + 2d^){ip), (17.15) 



which is the dual to the mysterious system ( |17.1| ) . As Wilson shows [28] , the system ( |17.14 ) 



and the corresponding hierarchy are bi-Hamiltonian and they possess a Miura map which is 
canonical. Taking the Symbol of the corresponding modified systems, we get the hierarchy 
of the dual-to-mysterious mKdV extensions. For N > 2, the corresponding Miura map 
comes out of the factorization. 

£, = (d-vi)---{d- VN)id + VN)---{d + vi). (17.16) 

Having thus realized the mysterious extension as a Symbol of a coupled 2-component 
KdV system, one may hope to find the remaining unresolved extension, the strange one, 
among known 2-component integrable systems extending the KdV. This was suggested to 
me by G. Wilson, and his advice was right on the mark. 

There are two known integrable super-extensions of the KdV hierarchy [22]: the super- 
KdV system [9] 

Ut = d {3u^ + Uxx + 3^^x) , (17.17a) 

6 = 3n^^ + 6uCx + 4.Cxxx, (17.17b) 



Dark Equations 



435 



and the supersymnietric KdV equation of Manin and Radul [21]: 
ut = d (3n^ + u^^ + 3^^x) , 

St ^ "JUj;^ + oUl^x + <,xxxi 

where ^ is now an odd variable ( "fermion" ) . 

Taking the Symbol of system (|17.18| ) we arrive at the strange system ( |6.29 ): 



(17.18a) 
(17.18b) 



ut = d (3u2 + Uxx), ^t = {d (3n + S^) ) {ip) . 
The fact that (p is a descendant of an odd variable is no longer material. 



(17.19) 



The system ( 17.18 ) is a Lax system, with the Lax operator [21] 



L = V* + TV, 

T = i-eu, V 



d_ 

de 



+ ed, 



0. 



(17.20a) 
(17.20b) 



This provides us, upon extracting the Symbol, with the whole strangely-extended KdV 
hierarchy. Also, the system ( |17.18 ) is Hamiltonian, as is the corresponding hierarchy, and 
they possess a Miura map that is canonical [23]; the strangely-extended mKdV hierarchy 
results thereby. (The strange system ( |17.19| ) for bosonic ip appears as an universal fluid- 
dynamical system, and also as a continuous limit of certain discrete mechanical spring 
models [24]. In the paper [24], Miller and Clarke consider the system 



ut 



a(3n^-Fu^^), V3t = (5(6A;m-F(9^))((^), 



(17.21) 



and conclude after a detailed analytical study that only for A; = 1 or A; = 1/2 does this 
system exhibit behavior that can be reasonably described as regular.) 

If we now extract the Symbol from the super-KdV system ( 17.17| ), we arrive at the 
Lax-type extension ([6.221): 



nt 



d (3n^ + Uxx) , 



n 



[?,ud + ?,du + Ad^) {if). 



(17.22) 



Since the super-KdV system ( 17.171) and the associated hierarchy are bi-Hamiltonian, with 
a canonical Miura map [9], we obtain the corresponding extension of the mKdV hierarchy; 
this also explains the unexpected at the time formula ( |8.28| ). 

The three examples above show that with the single exception of the linearized KdV 
extension, every other nontrivial KdV extension of Section 6 is either a Symbol of a nonlin- 
ear coupled integrable KdV system — or else its dual is. This suggests that our subject — 
linear extensions — in some large parts should be considered as a first approximation to 
the much more complex, difficult, and amorphous problem of nonlinear extensions. 

To test this point of view, let us examine the additional KdV extension (|9.9|): 



ut 



d (3n^ + Uxx) , 'ft = {aux + 2^9) (93). 



For a = 2, we get 

Ut = d {2,V? + Uxx) , 'ft = 2{u'p)a 



(17.23) 



(17.24) 
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and this is the Symbol of the Ito system [6] 

ut = d {3u^ + Uxx + V^) , '^t = '^iuf)x- (17.25) 

This system is bi-Hamiltonian [6]; the 2"'^ Hamiltonian structure is given by the Hamilto- 
nian matrix ( |9.23| ) for a = 2, and the f** Hamiltonian structure is 

d 
d 

Since the Miura ma p (|9.38D is no longer Hamiltonian, the existence of nonlinear mKdV 
extensions of which ( |9.41| ) is the Symbol (and which map into the Ito hierarchy), remains 
an open problem. 

It would be interesting to find nonlinear analogs of other linear extensions constructed 
in this paper, in particular for the Burgers, long waves, Benney, and KP hierarchies. It 
wouldn't be surprising if different nonlinear extensions have the same Symbol, as integrable 
systems are often realized as "invariant submanifolds" of many different larger ones. For 
example, the Burgers hierarchy results upon setting {h = 0} in the Dispersive Water Waves 
(DWW) hierarchy, mentioned in Section 12, while that same DWW hierarchy results in 
the KdV hierarchy upon setting {u = 0} for every second flow [12]. 

Remark 17.26. The search for nonlinear extensions can be made systematic. For the 
Hamiltonian case, most likely one has to take the Hamiltonian structure of linear extension 
and deform the Hamiltonian function itself. For the general case, one has to use the ansaltz 



ut = Y,<^s{v)+0[<f'^'), (17.27a) 

i+l 

n = Y, As{^) + O ((^^+2^ , £ G Z+ (17.27b) 

where Os{'-p) and As{'-p) are s-linear in c/?; for each such fixed £, one can define the com- 
mutativty of the flows modulo appropriate 0-terms. The homogeneuity requirement then 
guarantees that the sums on s in ( |17.27| ) are finite, so the whole problem can be auto- 
mated. (For systems of the type considered in Section 9, with the lowest order extensions 
of the form 

ut = X^{u), (17.28a) 

^t = A{u), (17.28b) 



the analog of an £*^ order extension (17.27) is 



ut 



Y,Os{v) + oU+A, (17.29a) 



=0 



(^i = ^A(<^)+0((^^+ij, ^^^^^ (17.29b) 

s=0 
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As a somewhat funny example, consider the dual to the linearized KdV extension 

ut = d (3n2 + u^^) , (ft = {6ud + d^) (93). (17.30a) 

It can be completed to the fully nonlinear system 

ut = d (3n^ + Uxx) , ^t = {Qud + d^) (ip) + 3ipx^xx + ^l, (17.30b) 

which appears as the addition to the KdV equation of another copy of it, but in the 
variable 

if = log{u). (17.30c) 

It's very likely that this particular doubling device (which can be made unique by the re- 
quirement that (ft is a differential polynomial in f and a (nondifferential) rational function 
in u which vanishes together with (p) applies to the whole KdV hierarchy. 

Remark 17.31. The Ito example (|17.24 ) suggests that "nonlinearization" may be possi- 
ble, if at all, only for some special values of the parameters of a linear extension. Formulae 
(17.28) and ( |17.30 ) offer an opportunity to construct and evaluate some sort of obstruc- 
tions to nonlinearization; we won't pursue this route here. Instead, here is another instance 
of such quantization. Consider the extended Hamiltonian matrix (12.29|) for a = 0: 





u 


h 


V 


u/ 





d 


-ipx/h 


h 


d 









B= h\ d I (17.32) 

95 \^x/h 

The Hamiltonian h is still a Casimir of the extended matrix, but the Hamiltonian u is no 
longer so; instead, there is a new Casimir 

H = hp^. (17.33) 

In addition, the matrix B ( |17.32| ) admits the following Hamiltonian automorphism: 

U = u + Xip, h = h, if = if, A = const, (17.34) 

a 99-extension of the identity automorphism of the unextended Hamiltonian matrix ( p.2.2D . 
This implies that the old linear extended hierachy in [/, /i, ip variables becomes new fully 
non/inear extended hierarchy in n, h, ip variables. For example, the original 1** flow ( |12.1 ), 
( 12.15 ) now takes the form 

ut = uux + hx + Xux<p, ht = [h{u + X'p)]x, (17.35a) 

ipt = {u + Xip)ipx. (17.35b) 

This is so because the old Hamiltonian {h?u + /i^) /2 turns into 

H2 = = h Xfiuip + I — hip IS a Casimir 1 . (17.3dj 



All this is only for the value a = 0. Whether other values of a allow nonlinearization, 
I don't know. Moreover, the above picture can be generalized into the moments space: 
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Proposition 17.37. Let B be the Hamiltonian matrix ( 13.1^ ) for a = 0; 



E 

B^^ = 0. 

(i) The map 

n 
fc=0 



nAn+m-id + dmAn+m~i, n,meZ+, 
(frcmAm-l/Ao, Bmp = -(fxnAn-i/Ao, 



Ak{X^) 



n—k 



n G Z4.; y^ = ^, 



(17.38a) 
(17.38b) 
(17.38c) 



(17.39) 



is a Hamiltonian automorphism of the Hamiltonian matrix B {17.3t); 

(a) In addition to the old Casimir H = Aq, the extended Hamiltonian matrix B ( 1 7. 3^ ) 
has also a new one; 



H = Ao^^ 



(17.40a) 



(Hi) The automorphism ( 17. 3^ ) covers the automorphism ( 17.34 ) with respect to the 
Hamiltonian embedding s An = hW^^An = hu^: 



J\.ri 



hU'\ A 



hu'- 



U = u + Xip. 



(17.40b) 



Proof, (i) A lengthy but straightforward verification, 

(ii) From formula ( 17.39| ) for n = 1, we have U = Ai/h = {Ai + A(/?^o)/^o = u + Xtp, 
and then 



E(:)^.(vr-'-E(l!) 



hu^{Xip) 



n—k 



h{u + XifY 



hW 



J\n. 



As a Corollary, the old linearly extended Benney system ( |13.14 ), (13.18): 

A.n,t = An+l,x + nAn-i, Ao.x, n G Z+, 

ipt = ifxAi/Ao, 
generates a fully nonlinear extension 

An,t = An+l,x + nAn-lAo^x + XipAn,x 

+ Xipx[{n + l)An-nAn-iAi/AQ], n G Z+, 
Lpt = ifxiAi + XipAo)/Ao; 

this is so because the old Hamiltonian H2 = {A2 + ^q) /2 turns into 

{A2 + Al) /2 = {A2 + Al) /2 + XAnp + {X^Aqlp^/2 is a Casimir) . 



(17.41a) 
(17.41b) 



(17.42a) 
(17.42b) 



(17.43) 



Naturally, when Uy = 0, so that An = f^ u^dy becomes An = hu^ and An = /q {u+Xip)^dy 
becomes An = h{u + Xtp)"', equations ( |17.42|) reduce to the system ( 17.35|) . 

Moreover still, the preceding picture can be generalized for the fully dispersive case: 
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Proposition 17.44. Let B be the extended KP Hamiltonian matrix ( 14-14 ) for the a = 0; 



K - V I" /^ " ^ /4 , rT 

Km = j^v^. i:(-^) ys + l 



i-dri ^ u. 



n+m—fj. 






m— 1— S) 



E 



n \ 1 - 



(i) The map 

An = ^i ^ j AkQn-k{>>'p), nGZ+; (^ = (^ 



, n, ?TT.sZ_|_, (17.45a) 
(17.45b) 
(17.45c) 

(17.46) 



A:=0 



is a Hamiltonian automorphism of the Hamiltonian matrix B ( 17. 4t ); 

(a) In addition to the old Casimir H = Aq, the extended Hamiltonian matrix B ( 17.4^ ) 
has also a new one, 



H = Ao^^; 



(17.47) 



(Hi) The automorphism ( 17. 4t ) covers the automorphism ( 1 7. 34 ) with respect to the 
Hamiltonian embeddings 



■An = hQn{U), An = hQn{u) : 

An = hQn{K), An = hQn{u) =^ U = U + Xif. 



(17.48) 



Proof, (i) A lengthy but straightforward verification, requiring at the last step the identity 
([12, p. 66]) 






(17.49) 



(ii) Is obvious; 

(iii) Is based on the addition formula (a differential analog of the Newton binomial) 



Qniu + v) = J2[l)Q 



k{u)Qn-k{v), n e Z+. 



(17.50) 



which is easily proved by induction on n. 

As a Corollary, the old linearly extended KP system (for the flow Lt = [P+ , L] , P 



L'/2,L = d+Y.d-'-^Ai): 

j=0 



A .-^W _i4(2) + Vr "^ \a ^ 4^^+^) uG 

^n,t — ^n+1 2" ^l U+1 / ^"--^-/J-^O ' "- ^ 



V^i = V'x(Ai-4'V2)Mo, 



(17.51a) 
(17.51b) 
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becomes the fully nonlinear system 



A^,t = A« 1 - ^4'^ + E ( ;, + 1 ) ^-1-/^ (^0 + A^./2)(''+^) + A(A„^) 



(1) 



+ A E ( ll)\ An-^'f^^^^^ - ^n-l-M [p.A^ - Af (2)Mo 

(ft = (f^ (^Ai - a'^^/I + A(^^o) Mo; 
this is so because the old Hamiltonian (^2 + ^0) /^ turns into 

{A2 + Al) /2 = [A2 + 2AiXip + ^0 (AV^ + A(/^x) + vlg] /2 

= (^2 + ^0) /2 + A(Ai(/3 + ^o93x/2) + {X^Ao(p'^/2 is a Casimir) 



M 



(17.52a) 
(17.52b) 



(17.53) 



18 Scalar extensions associated with scalar Lax operators 
of order > 4 



Let -L be a scalar Lax operator of order A^ > 4: 

N-2 



j=0 
ind n„j +!-,„ ord 



Consider the 2 and the 3'' flow of the hierarchy 
(Li=) X„.iL) = [P+,L], P=(l'/''" 



j2/N 



92 + |rUo9 + 



3 /3-7V 



-uo,x + ni 



N N \ 2 

(The last expression results from the condition 

ord([P+,L]) <N -2. 
The motion equations (|18.2D , at least those we shall need below, are: 

^2 (no) = no,i = (2 - N)u^^^ + 2u^^\ 

, X (iV-l)(7V-2) (3) (2) fi) 2(7V-2) 



(1) 



-^^3(^*0) = ^io,t 



7V-3V^f3) 3(A^-3)^(2) , , (1) 3(iV-3)^, ^(1) 



uf ' + 3u^ 



2 ; '^ 2 

Let's look for scalar extensions of the flows X2 and X3: 
(Xf*) ipt = A2{^) = (ano + 79') (<^), 
(A3°^*) (^t = Asi^) = ((a4^^ + 6ni) + /no5 + 59=^)('^), 



AT 



-UqUq 



(18.1) 

(18.2) 
(18.3) 
(18.4) 

(18.5) 

(18.6a) 
(18.6b) 

(18.7) 

(18.8a) 
(18.8b) 
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where a, 7; a, b, /, g are unknown constants. The commutativity condition ( |1.13 ) 

X2 (au\) ' + bui + fuodj - X3{auo) = [A2, A3], 
yields the relations: 



(18.9) 



27/ = Sag, 

'3{N -3) 



a 



N 



+ f 



2{N -2) 
' N -■ 



jf-3ag + 2^a = f{2-N), 

26 = 3a, 

..s ,(^-l)(^^-2) /iV-3 

-ag + ja = a(2 - N) - b^ ^-^ ^ - a 



76 = 2a + 6 + 3a- 



N -3 



If 



a = => a = b = f = 
is the decomposed case. So let us consider 
a/0. 



(18.10a) 

(18.10b) 

(18.10c) 
(18.10d) 
(18.10e) 

(18.10f) 
(18.10g) 

(18.11) 
(18.12) 



If 7 = then (|18.10a| ) yields 5 = 0, ( |18.10dD yields / = 0, and |( |18.10bD & ( |18.10e| )} 
force a contradiction. 
So, let 



a/0, 7/0. 

Then ( |18.10e| ) yields 
, 3a 



and from ( 18.10d ) we get 

307 
■' 2 

Substituting this into ( 18.10a| ), we find 

9 = 1^- 
Then ( |18.10c| ) yields 



(18.13) 



(18.14) 



(18.15) 



(18.16) 



a = ^(2 + 7-iV), 



(18.17) 
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while ( |l8.10b| ) provides 
2 



a 



jN' 



Combining ( 18.18| ) with the preceding formulae ( 18.14D -( |18.17l ), we find 

3 



27 2 . 3 37 

"=iV' ^ = ^' ^=iV' ^=iV' "=2^ 



(2 + 7 -TV). 



(18.18) 



(18.19) 



Substituting this into (18.10^), we get 



or 



Thus, 



,37 3 , ,,, 3-2 Af- 3 



(7 - 1)72 = 2 + 7-iV + iV-3 = 7-l 
^ (7-1) (72 - 1) = 0. 

7 = ±1. 



For 7 = 1, formulae (18.19|) yield 



2 ^3,3 3 3-iV 



(18.20) 
(18.21) 

(18.22) 



and formulae ( |18.3| ), (|18.4| ) show that this is precisely the Lax-type extension. The case 
7 = — 1 is the dual extension. 

Thus, scalar Lax operators of order > 4 have just two 1-component extensions: the Lax- 
type one and its dual. Extensions with 2i components, i £ N, come from Lax equations 
with Lax operators of the form 



N-2 



L-t = d^+Y, n,d^-^-' + Y, Vsd-' 



^i+s 



(18.23) 



1=1 



s=l 



(see § 3.5 in [16]). If L is self-adjoint then we can take 



L-t = L + ^65-i6, 



where the ^'s are odd. If L is skew-adjoint, we can take 



(18.24) 



L-t = L + ^V'.5-V. 



s=l 



where the 93's are even. The simplest case of the latter sort is 
L = d^ + ud + du, 



(18.25) 



(18.26) 
(18.27) 
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19 Composition of Hamiltonian extensions 

Suppose we have a Hamiltonian matrix B = B{q} = {Bnm) in g-variables, and suppose 
we have two Hamiltonian extensions of B, B^ and B^, of the form 



B' 



B' 





qm 


"PP 


qn 


-Dnm 




fa 


\ ^am 


1 




qm 


i^u 


Qn 


( Bnm 




i^/i 


\ ^^m, 






il9.l) 



(19.2) 



where the operators Oam depend only on {(^}, and operators O^m depend only on {ip}- 
Proposition 19.3. If B^ and B^ are Hamiltonian then so is their direct sum 



B 





qm 


fp 


Viy 


( n 

q-n -Djim, 






fa 


Oam 








^M 


\0^m 





/ 



(19.4) 



Proof. By the main result of the Hamiltonian formalism [14, p. 49], to show that a matrix 
is Hamiltonian we need only verify the Jacobi identity for arbitrary linear Hamiltonians: 

{H,{F,G}} + {F,{G,H}} + {G,{H,F}}^^, y H,G,F. (19.5) 

1) If all Hamiltonians H, F, G, are g-independent, {H, {F, G}} = 0, V i7, F, G; 

2) If two out of three Hamiltonians H, F, G are ^-independent, {H, {F, G}} = 0, 
V H, F, G; 



3) If all Hamiltonians H, F, G are (p- and ^'-independent, then (|19.5| ) is the criterion 
for the extendee matrix B to be Hamiltonian; 

4) It leaves us with the case when two out of three Hamiltonians are cp- and V-inde- 
pendent, and the third one is g- independent. If that third Hamiltonian is -(/'-independent, 
of the form 



f/V, 



(19.6a) 



then (19.5) is the condition for the extended matrix B^ to be Hamiltonian; if that third 
Hamiltonian is 99- independent, of the form 



V'lP, 



(19.6b) 



then ( 19.5 ) is the condition for the extended matrix B^ to be Hamiltonian; in general, any 
linear g- independent Hamiltonian is a sum of two terms, of the type ( |19.6 ), and the LHS 
of ( |19^ is trilinear in H, F, G. ■ 
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